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ABSTRACT 

Estimates  of  the  sampling  frequency  required  to  describe  armual  and  seasonal 
(summer  stratification  and  ice-covered  period)  chemical  concentrations  in  lake  water  within 
a  given  tolerance  of  the  true  mean,  with  95%  confidence,  were  made  using  data  collected 
for  three  lakes.  Visual  inspection  and  Monte  Carlo  simulation  of  the  data  suggest  that 
most  parameters  approximate  a  normal  distribution  and  thus  simple  parametric  statistical 
methods  may  be  applied  to  lake  data.  Graphs  are  presented  showing  95%  confidence 
intervals  from  seasonal  means  for  sample  sizes  of  from  one  to  eight  for  12  chemical 
parameters.  These  estimates  will  be  useful  in  planning  lake  water  quality  surveys.  The 
effect  of  sampling  season  is  also  discussed. 


INTRODUCTION 

In  order  to  assess  the  extent  of  the  aquatic  resources  sensitive  to  acidification,  it  is 
necessary  to  characterize  lakes  over  a  wide  geographic  range  with  respect  to  important 
parameters  such  as  alkaUnity  (acid  neutralizing  capacity),  pH,  SO^,  and  base  cations 
(Harvey  et  al.  1981,  Dillon  etai.  1984,  Haines  1981).  Characterization  of  lakes  at  points 
in  time  is  also  essential  information  needed  to  assess  temporal  changes  in  the  chemistry  of 
lakes.  It  is  particularly  important  to  establish  baseline  chemical  conditions  of  a  large 
number  of  lakes  in  order  that  future  changes  can  be  readily  evaluated.  However,  to  place 
statistically-defined  confidence  about  any  chemical  characterization,  it  is  necessary  to  collect 
more  than  one  sample.  On  the  other  hand,  it  is  obviously  desirable  to  collect  as  few  as 
possible  samples  with  minimal  reduction  in  information.  The  purpose  of  this  paper  is  to 
define  the  minimum  sampUng  frequency  necessary  to  adequately  define  lake  chemistry  with 
respect  to  the  acidification  phenomenon. 

This  subject  has  been  discussed  in  general  in  many  papers  (eg.  Ellis  1980.  Ellis  and 
Lacey  1980,  Walker  1982,  Cotter  1985)  and  texts,  but  the  specific  case  of  acidification  of 
lakes  has  not  previously  been  addressed.  To  estimate  the  required  sample  size,  past 
estimates  of  the  variability  of  a  given  water  body  are  required.  In  addiuon,  most  statistical 
methods  recommended  for  water  quality  surveys  require  the  data  fit  a  normal  distribution. 
Estimation  of  variance  and  normality  requires  relatively  large  numbers  of  data.  In  many 
cases,  acquisition  of  this  information  is  impractical  or  not  economically  feasible. 


In  this  paper,  we  utilize  an  existing  data  base  produced  as  a  result  of  an  intensive 
sampling  effort  on  three  representative  lakes  to  develop  the  statistical  considerations 
necessary  for  planning  a  survey  of  similar  lakes  with  respect  to  sensitivity  to  acidification. 
This  included  1)  using  the  estimated  variability  of  the  annual  chemistry  of  the  three 
representative  lakes  to  determine  the  minimum  sampling  frequency  necessary  to  adequately 
define  lake  chemistry  and  2)  assessing  the  relative  accuracy  of  1  random  sample  to  estimate 
the  annual  mean.  In  addition  the  data  were  used  to  3)  estimate  the  distribution  of  the 
selected  water  parameters  and  test,  using  a  Monte  Carlo  simulation,  the  normal  distribution 
assumption  and  thus  the  applicabihty  of  parametric  statistical  methods.  Since  many 
sampling  programmes  are  designed  to  estabhsh  aimual  means  but  are  restricted  to  sampling 
only  seasonally,  we  used  our  data  set  to  4)  test  whether  seasonal  means  were  significantly 
different  from  annual  means,  and  5)  assess  the  potential  bias  with  estimating  annual  means 
from  one  seasonal  sample. 


STUDY  SITE 

The  lakes  discussed  here  have  been  intensively  studied  since  1976.  All  are  in  the 
Muskoka/Haliburton  region  of  Ontario,  about  200  km  north  of  Toronto  and  200  km 
southeast  of  Sudbury,  the  closest  centres  of  industrial  activity.  They  are  in  a  zone  of  high 
deposition  of  atmospherically-transported  strong  acid  and  sulphate  (Dillon  et  al.  1978, 
Likens  and  Butler  1981,  Dillon  etal.  1988),  but  none  receive  pollutants  directly  from  point 
sources.  The  three  lakes  discussed  here,  Blue  Chalk,  Chub  and  Harp,  are  relatively  small 
(Table  1),  but  typical  in  size  of  the  approximately  12,000  lakes  in  this  region  (Cox  1978). 


Their  use  is  primarily  recreational  with  seasonally-occupied  dwellings  on  the  shores  of  all 
3  lakes. 

The  study  lakes  are  situated  on  the  Precambrian  Canadian  Shield,  in  a 
heterogeneous  area  of  metamorphic,  plutonic  and  volcanic  silicate  bedrock.  The 
unconsolidated  surficial  deposits  in  these  lakes 'watersheds  consist  of  sub-glacially  deposited 
tills  of  variable  thickness,  with  some  accumulation  of  sand  and  gravel  deposited  by  glacial 
meltwaters,  and  of  organic  matter  (peat)  in  poorly  drained  areas.  The  catchments  of  the 
3  lakes  each  have  areas  with  relatively  thick  till,  with  deposits  in  the  Harp  Lake  catchment 
known  to  reach  >  10m.  Thus,  the  3  lakes  are  considered  to  be  somewhat  less  sensitive  to 
acidification  thzin  those  lakes  in  the  high  deposition  zone  with  entire  catchments  covered 
by  only  very  thin  till. 

Lakes  similar  to  these  three  (or  less  sensitive)  constitute  -70%  of  the  sample  of 
about  2000  lakes  surveyed  in  the  high  deposition  portion  of  Precambrian  southern  Ontario 
(Neary  et  al.  1989).  Lakes  that  are  not  highly  sensitive  to  acidic  deposition  were 
deliberately  chosen  for  this  analysis  to  minimize  or  eliminate  long-term  trends  in  chemistry. 
To  test  this,  time  trend  analysis  was  used  (reported  elsewhere:  Dillon  el  al.  1987,  Dillon 
unpubl.  studies)  to  demonstrate  that  there  were  either  no  or  extremely  slight  (in  the  case 
of  sol')  long-term  trends  in  the  chemical  parameters  reported  in  this  study  for  any  of  the 
3  lakes  over  the  years  in  question. 


METHODS 

i)          Parameter  selection 

In  order  to  characterize  lakes  with  respect  to  their  sensitivity  to  acidification,  12 
parameters  were  considered  important. 

Alkalinity,  the  acid  neutralizing  capacity  (ANC)  of  the  water,  is  probably  most 
important,  particularly  since  it  is,  unlike  pH,  conservative  with  respect  to  CO2 
concentrations.  Biological  organisms,  however,  respond  to  ambient  pH  rather  than  ANC 
(Mierle  et  al.  1986).  The  SO*  content  is  an  indication  of  the  importance  of  long-range 
transport  of  sulphuric  acid  and  its  precursors,  since  the  'background'  or  natural  SO4  level 
in  lakes  in  this  region  is  low  (Neary  and  Dillon  1988).  It  is  also  a  mobile  anion  with 
respect  to  flux  through  lakes'  catchments,  and  this  controls  to  a  large  degree  the  extent  and 
rate  of  acidification  of  aquatic  systems  (Reuss  and  Johnson  1986).  The  cations  Ca,  Mg,  Na, 
and  K  are  products  of  weathering  and  ion  exchange  processes,  which  are  believed  to  occur 
at  an  increased  rate  in  areas  affected  by  acidic  deposition.  Dissolved  organic  carbon  (DOC) 
and  colour  aie  strongly  related  to  the  organic  acid  content  of  natural  waters,  and  indicate 
qualitatively  the  importance  of  naturally  occurring  organic  acids.  Aluminum  is  important 
because  of  its  toxicological  effects  on  fish  and  other  organisms.  Chloride  is  often  chosen 
as  a  conservative  element  generally  believed  to  be  largely  imaffected  by  acidic  deposition, 
and  is  only  rarely  a  significant  component  of  the  atmospheric  acids.  Conductivity,  a 
measure  of  the  total  ionic  strength  of  the  water,  provides  a  general  measure  of  how  'dilute' 
the  water   is.      Although   nitrate   and   ammonium   are   significant   components  of  the 


atmospheric  deposition  in  eastern  North  America,  they  are  not  included  here  because  of 
their  biological  reactivity;  levels  measured  in  lakes  are  not,  in  general,  indicative  of  input 
rates. 

ii)         Sampling  Frequency 

Water  samples  were  taken  for  chemical  analysis  from  Blue  Chalk,  Harp  and  Chub 
Lakes  from  1976  to  1985  as  part  of  an  investigation  of  these  three  lakes  and  their 
watersheds.  The  lakes  were  generally  sampled  bi-weekly  (twice  per  month),  more 
frequently  at  spring  and  fall  overturns,  and  less  frequently  in  winter,  with  the  exception  that 
during  the  years  1977  to  1981,  the  lakes  were  visited  every  week  during  the  summer  period. 
The  number  of  samples  taken  under  ice  cover  varied  from  3  to  8  with  4  visits  being  the 
average,  Le.  approximately  one  sample  per  month.  The  lowest  number  of  samples  per  year 
was  20  and  21  visits  per  lake  in  the  hydrologic  years  1982-83  and  1983-84.  The  most 
frequently  sampled  hydrologic  year  was  1977-78  with  31,  35  and  36  visits  respectively  to 
Blue  Chalk,  Harp  and  Chub  Lakes. 

iii)        Sampling  Technique 

Samples  were  collected  at  the  deepest  spot  in  each  lake.  A  peristaltic  pump  and 
Tygon  tubing  were  used  to  obtain  water  at  2  m  intervals  commencing  at  1  m  below  the 
surface  to  1  m  above  the  bottom.  Aliquots  from  all  depths  were  combined  such  that  a 
single  volume-weighted  composite  sample  was  obtained  for  each  of  the  epilimnion, 
metalimnion  and  hypolimnion  during  the  period  of  thermal  stratification,    lliese  three 


samples  were  analyzed  separately  for  chemical  parameters.  Under  conditions  of  thermal 
homogeneity,  water  from  all  depths  was  pooled  to  give  a  single  volume-weighted  sample 
representative  of  the  entire  lake.  Sample  containers  were  filled  by  pouring  the  composite 
water  through  80  uM  Nitex  mesh  to  remove  coarse  particulate  material.  A  complete 
description  of  sampling  methods  is  available  in  Scheider  et  al.  (1983)  and  Locke  and  Scott 
(1986). 

iv)        Analytical  Methodology 

Chemical  analyses  were  performed  at  Dorset  or  Toronto,  depending  on  the 
perishability  of  the  parameter.  All  analytical  techniques  are  described  in  Ontario  Ministry 
of  the  Enviroimient  (1983)  and  the  laboratory  performance  is  reported  in  Ontario  Ministry 
of  the  Environment  (1986).  The  analytical  method  used  and  the  average  analytical  standard 
deviation  based  on  replicate  samples  and  standards  are  summarized  in  Table  2  for  the  12 
chemical  parameters  considered  in  this  paper. 

The  methodology  for  7  of  the  12  parameters  (Al,  Ca,  Mg,  Na,  K,  Conductivity,  pH) 
did  not  change  significantly  over  the  course  of  this  study.  Method  alterations  for  3 
parameters  (CI,  DOC,  and  Colour)  did  not  result  in  any  systematic  differences  before  and 
after  the  changes,  as  shown  by  intercomparison  of  methods.  Alkalinity  was  measured  by 
titration  to  the  equivalence  point  (Gran  titration)  from  July  1979.  Prior  to  this  date, 
titration  to  a  fixed  end-point  (pH  4.5)  that  was  lower  than  the  equivalence  point  (usually 
pH  4.8-5.1)  was  carried  out  routinely.  Although  there  is  a  defined  relationship  between  the 


two  titration  results  (Harvey  etal.  1981,  Henriksen  1982),  we  have  used  only  those  data 
collected  by  the  Gran  Method. 

Sulphate  was  measured  by  the  methylthymol  blue  (MTB)  method  to  June  1980  and 
by  ion  chromatography  thereafter.  Because  of  a  significant  colour  interference  mth  the 
MTB  method  (Kerekes  1984),  only  the  results  obtained  after  June  1980  are  considered 
here. 

v)        Statistical  Methods 

Data  Reduction 

The  analytical  results  of  the  individual  samples  collected  on  each  date  were  edited 
to  delete  obviously  incorrect  values  using  the  principles  of  charge  balance,  the  relationship 
between  ionic  strength  and  conductivity  and  the  theoretical  relationships  between  pH, 
alkalinity  and  DIC  (Stumm  and  Morgan  1981).  DIC  was  measured  on  all  samples  but  the 
results  are  not  reported  here  (P.  Dillon,  unpubl.  studies).  Less  than  2%  of  the  results  were 
discarded  in  this  fashion. 

During  periods  of  lake  stratification,  the  chemical  results  from  each  strata  were 
volume-weighted  according  to  lake  morphometry  to  give  a  single  value.  This  data  was  then 
pooled  annually  (AN)  and  over  the  whole  ten  year  study  (TOT).  For  some  statistical 
analy.ses,  the  data  were  also  pooled  by  functionally-defined  seasons  over  the  whole  ion- 


year  study,  including  spring  overturn  (SO),  summer  stratification  (SS),  fall  overturn  (FO), 
ice  free  (IF  -  which  includes  all  of  the  above)  and  ice  covered  (IC). 

Sample  Distribution 

To  determine  the  number  of  samples  required  for  a  given  level  of  statistical 
confidence,  the  underlying  variance  of  the  sampled  population  has  to  be  estimated.  As 
there  is  significant  year-to-year  variability  in  these  data  sets  (see  below),  the  best  overall 
estimate  of  variance  about  an  annual  (or  seasonal)  mean  would  be  the  average,  over  ten 
years,  of  the  individual  annual  (or  seasonal)  variances.  However,  we  decided  to  be 
conservative  and  to  over-estimate  the  variance  by  using  the  variance  obtained  from  the  ten- 
year  pooled  data  sets  (TOT,  SO,  SS,  FO,  IF  and  IC). 

None  of  these  samples  were  collected  randomly  over  time.  Rather  a  systematic 
sampling  strategy  was  employed,  primarily  because  of  the  expected  seasonal  trends  in  some 
parameters.  These  seasonal  trends  as  well  as  the  significant  within-year  similarity  (as 
opposed  to  between  years)  causes  serial  dependence  in  the  data.  When  all  samples  are 
lumped  together,  dependencies  of  this  sort  can  distort  the  sampling  distribution  and  make 
the  estimation  of  sample  size  for  a  given  confidence  level  problematic.  Nonetheless,  we  will 
assume  that  the  ten-year  pooled  data  sets  consist  of  approximately  independent  and 
identically  distributed  variates. 

Given  the  variance  of  a  Normal  sample  distribution,  the  calculation  of  required 
sample  size  is  straightforward.     However  if  the  distribution  is  non-normal,  alternative 


procedures  must  be  used.  We  tested  the  fit  of  the  annual  data  and  the  two  non-overlapping 
seasonal  data  sets  with  the  greatest  sample  size  (SS  and  IC)  to  the  normal  distribution  using 
Lilliefor's  test  (Conover  1971).  In  addition,  we  compared,  for  a  given  sample  size,  the 
estimated  confidence  intervals  about  the  annual  mean  obtained  with  the  normal 
approximation  to  those  obtained  with  a  Monte  Carlo  approximation  that  makes  no 
assumption  about  sample  distribution. 

The  Normal  Approximation 

The  normal  approximation  for  an  interval   'L'  about  the   mean,  within  which 
observations  will  fall  95%  of  the  time  (Snedecor  and  Cochran  1980)  is: 

r 

L  =  1.96  S  /  v/n 

where  S  =  population  standard  error 
n  =  number  of  samples 

We  used  the  sample  standard  deviations  as  an  approximation  of  the  population  value.  The 
simple  calculation  for  L  was  made  on  n  =  1  to  8  for  the  aimual  data  (AN)  and  the  two 
seasons  (SS  and  IC). 

For  comparison  a  binomial  test,  which  makes  no  assumptions  about  the  distribution 
of  the  samples,  was   used  to  estimate  sample  size  (Cotter  1985Siumm 


Monte  Carlo  Comparison 

The  95%  confidence  interval  'L'  estimated  with  the  normal  approximation  was 
compared  with  a  Monte  Carlo  simulation.  The  simulation  drew  every  possible  combination, 
without  replacement,  of  the  specific  sample  size  (n  =  1  to  3).  The  distance  'L'of  the  95% 
confidence  interval  was  set  as  an  upper  and  lower  boundary  about  the  mean  and  the 
number  of  simulated  means  falling  beyond  each  boundary  was  counted.  Expressed  as  a 
percent  of  the  total  number  of  cases,  the  count  was  compared  to  the  normal  95% 
confidence  interval. 


RESULTS 

i)         Descriptive  Chemistry 

1)        Annual  Means 

The  three  lakes  represent  a  range  of  chemical  composition  normally  observed  on 
Precambrian  terrain  in  much  of  Ontario  (Table  3).  All  are  dilute  lakes  (mean  conductivity 
over  10  years  ranged  from  29.7  to  35.3  ^S)  with  low  levels  of  major  ions.  Blue  Chalk  lake 
is  a  very  clear  lake  (colour  =  6.0  Hazen  units)  with  low  dissolved  organic  carbon  (DOC  = 
1.8  mg  r^)  and  moderate  alkahnity  (4.0  mg  1'^  as  CaCOj).  Harp  Lake  is  intermediate  in 
color  (21.1  Hazen  units)  and  DOC  (4.1  mg  1'^)  and  also  low  alkalinity  (3.2  mg  1'^  as 
CaCOj).   Chub  lake  is  a  coloured  (dystrophic)  lake  (46.5  Hazen  units)  with  a  DOC  value 
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of  4.9  mg  r^  and  very  low  alkalinity  (0.70  mg  1'^  as  CaCOa).   Chub  Lake  has  the  highest 
average  Al  concentration  (96  /ig  1"^)  and  lowest  average  pH  value  (5.66)  of  the  study  lakes. 

2)         Seasonal  Variability 

The  average  seasonal  values  and  the  annual  average  for  each  parameter  for  the  three 
lakes  are  shown  in  Fig.  1-3.  The  measured  standard  deviations  and  means  for  each 
parameter  for  annual  and  each  of  two  seasons  (SS  and  IC)  are  summarized  in  Table  3.  In 
addition,  the  aimual  and  seasonal  means  for  each  year  were  matched  and  compared  over 
all  ten  years  using  the  Wilcoxon  matched  pairs  signed  rank  test  and  are  shown  in  Table  4. 

There  appears  to  be  little  seasonal  bias  compared  to  the  annual  means  as  only  a  few 
parameters  showed  any  significant  (p<0.05)  seasonal  trend  (Fig.  1  to  3  and  Table  4).  Most 
of  the  differences  between  the  long  term  seasonal  and  annual  means  were  less  than  the 
analytical  standard  deviation  (Table  4). 

The  summer-stratified  and  ice-cover  seasons  encompass  the  majority  of  the  sample 
dates  used  to  calculate  the  annual  mean  and,  not  surprisingly,  both  are  close  to  the  annual 
mean  value  (Table  3).  The  ice  free  season  (IF)  is  included  here  because  lake  surveys  may 
be  conducted  in  which  the  thermal  condition  in  specific  lakes  are  not  measured.  It  appears 
that  there  is  no  significant  difference  between  IF  and  the  annual  means  for  most  of  the 
parameters  in  the  three  lakes  (Table  4).  Some  general  trends  are  apparent,  however.  Na, 
and  possibly  Mg,  concentrations  tended  lo  be  greater  than  the  annual  mean  during  most 
SS  seasons  in  all  three  lakes  and  less  during  the  IC  season.   I'ive  ot  the  12  parameters  were 
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consistently  lower  or  higher  than  the  annual  mean  during  the  IC  season  in  highly  coloured 
Chub  Lake.   Some  specific  seasonal  trends  of  differences  will  be  discussed  further. 

For  conductivity  the  seasonal  differences  were  well  within  the  standard  deviation  (sd) 
of  the  analytical  methods.  All  three  lakes  had  their  lowest  values  during  the  spring  overturn 
(only  Blue  Chalk  significantly  different,  A  =  -0.7  /iS)  and  the  highest  values  during  ice  cover 
(only  Chub  lake  significant,  A  =  0.7  fiS). 

The  pH's  of  all  three  lakes  were  consistently  less  than  the  annual  average  under  ice 
cover  and  during  spring  overturn.  Blue  Chalk  lake  (pH  =  6.59)  has  consistently  lower 
values  during  IC  (A  =  -0.17)  and  consistently  higher  values  over  the  summer  season  (A  = 
0.04).  In  contrast,  Chub  Lake  (pH  =  5.66)  had  consistently  lower  values  during  fall 
overturn  compared  to  the  annual  mean  (A  =  0.09). 

Calcium,  CI,  Na,  K  and  Mg  seasonal  mean  values  were  within  the  analytical  standard 
deviation  of  the  annual  mean  value.  However,  despite  the  small  differences  in  means  there 
is  some  seasonal  trend  of  consistently  lower  Na  concentrations  during  SS  and  higher  Na  and 
Mg  concentrations  during  IC  compared  to  annual  means  over  the  study  period. 

For  sulphate,  the  difference  between  the  seasonal  means  and  annual  means  were 
generally  within  the  sd  of  both  analytical  methods  (Table  3).  There  is  a  bias  to  lower  than 
average  annual  values  at  spring  overturn  (A  =  -0.10  to  -0.15  mg  1"^)  and  the  highest  values 
of  the  year  occurred  during  IC  (a  =  0.102  to  0.270  mg  1"^). 
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DOC  varied  little  throughout  the  year  for  Blue  Chalk  and  Harp  Lakes,  remaining 
within  the  standard  deviation  of  the  analytical  sd  (0.148  mg  1"').  In  Chub  Lake  the  mean 
DOC  concentrations  were  consistently  lower  than  during  spring  (A  =  -0.33  mg  "^)  and 
consistently  higher  during  ice  cover  (A  =  0.24  mg  1'^)  compared  to  the  aimual  mean. 

Colour  varied  between  seasons  differently  in  the  three  lakes.  Although  there  was 
no  significant  trend,  low  coloured  Blue  Chalk  Lake  (6.0  Hazen  units)  had  lower  colour 
under  ice  cover  in  each  of  the  3  years  it  was  measured  and  somewhat  lower  values  during 
the  spring  overturn.  Harp  lake,  with  intermediate  colour  (21.1  Hazen  units)  had 
consistently  higher  color  during  IC  and  slightly  higher  values  during  FO.  Highly  coloured 
Chub  Lake  (46.5  Hazen  units)  had  the  lowest  values  during  SO  and  highest  values  during 
FO. 

ii)        Data  Distribution 

1)         Histograms 

Histograms  of  the  multi-year  data  set  are  presented  for  each  parameter  for  all  three 
lakes  (Fig.  4  to  6).  Most  of  the  parameters  deviated  from  normality  according  to  Lilliefor's 
test  (Table  5).  The  theoretical  normal  distribution  about  the  mean  is  marked  on  the 
histograms  for  comparison.  Visual  inspection  reveals  that  most  of  the  parameters  have 
distributions  which  approximate  the  theoretical  normal  distribution  or  tend  lo  be  leptokurtic 
(values  over-represented  around  the  mean).    There  are  no  apparent  trends  in  deviation 
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from  normality  for  any  one  parameter  or  within  any  one  lake.  Specific  cases  are  described 
below. 

Conductivity  measurements  had  a  leptokurtic  distribution  about  the  mean  in  Blue 
Chalk  and  Harp  Lakes.  Chub  Lake  has  a  bimodal  distribution  with  one  peak  at  27.5  /iS 
and  another  at  30  /iS. 

Alkalinity  measurements  in  Chub  Lake,  which  had  the  lowest  mean  (0.70  mg  CaCOj 
r^),  were  negatively  skewed,  le.  a  larger  number  of  values  less  than  the  mean  with  a  long 
tail  above  the  mean.  Blue  Chalk  and  Harp  Lakes,  with  moderately  low  alkalinity,  have 
values  distributed  more  closely  to  the  normal-shaped  curve. 

pH  measurements  were  somewhat  positively  skewed  for  Blue  Chalk  Lake,  while  they 
were  close  to  normal  or  slightly  leptokurtic  in  Harp  and  Chub  Lakes. 

Calcium  measurements  appear  to  be  close  to  normally  distributed  in  all  three  lakes. 
Chub  Lake  values  are  slightly  more  platykurtic  and  the  other  two  are  slightly  leptokurtic 
in  shape.  Magnesium  is  perhaps  the  most  evenly  distributed  parameter  and  the  distribution 
is  close  to  the  normal  curve  in  all  three  lakes.  Sodium  concentrations  tended  to  be 
negatively  skewed  in  Harp  Lake  but  leptokurtic  in  the  other  two  lakes.  Similarly,  potassium 
has  a  very  peaked  distribution  in  Blue  Chalk  and  Harp  Lakes.  The  stair-casing  of 
alternating  high  and  low  frequencies  over  the  range  of  measurements  for  the  above 
parameters  is  probably  an  artifact  of  the  roimding  method  used  in  reporting  the  values. 


14 


The  chloride  concentrations  are  close  to  normally  distributed.  Blue  Chalk  Lake  had 
the  most  peaked  distribution  with  the  lowest  variability. 

Sulphate  measurements  had  slightly  different  distributions  in  each  of  the  three  lakes. 
The  Chub  lake  data  appears  to  be  bimodal,  with  one  peak  at  9.5  and  the  other  around  7.5 
mg  r\  In  Harp  Lake  the  measurements  appeared  to  be  platykurtic  about  the  mean  while 
they  appeared  to  be  normally  distributed  in  Blue  Chalk  Lake. 

Harp  and  Blue  Chalk  Lakes  have  leptokunic  distributions  of  DOC  compared  to  the 
normal  curve,  with  a  slightly  negative  skew  in  Harp  Lake.  Chub  Lake  DOC  concentrations 
are  platykurtic  in  distribution. 

For  colour  measurements,  there  are  generally  more  values  at  the  mean  than  expected 
by  the  normal  distribution.  This  is  most  pronounced  in  Blue  Chalk  Lake,  which  has  the 
lowest  mean  colour  of  the  three  lakes. 

There  is  a  large  variability  in  aluminum  measurements  when  all  the  data  is 
considered.  Due  to  recent  improvements  in  analytical  methods,  data  collected  after  1984 
will  be  a  more  accurate  and  precise  measurement  of  water  concentration.  For  the  entire 
data  set,  Al  concentrations  in  Blue  Chalk  Lake  are  strongly  negatively  skewed,  while  .A.1 
concentrations  in  Harp  Lake  are  only  slightly  negatively  skewed.  The  distribution  of  Al 
values  is  leptokurtic  in  shape  in  Chub  Lake. 
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2)         Monte  Carlo  Comparison 

The  results  from  the  count  of  Monte  Carlo  simulation  cases  falling  beyond  the 
boundaries  of  the  estimated  tolerances,  with  95%  confidence,  for  a  given  sample  size  are 
shown  in  Tables  6a  to  61.  The  results  for  all  12  parameters  generally  have  counts  less  than 
5%,  i^e.  very  close  to  the  95%-quantile  falling  within  the  parametric  confidence  boundary. 
These  results  tend  to  verify  the  normality  assumptions  and  that  the  parametric  statistics  are 
robust  enough  to  be  applied  to  this  lake  chemistry  data.  The  highest  count  beyond  the  95% 
C.I.  in  the  Monte  Carlo  simulation  was  10.3%  (Na,  Harp  Lake,  IC,  Sample  size  n  =  1). 
The  lowest  counts  were  0%  of  the  cases  falling  beyond  the  95%  C.I.  This  occurred  in  a 
number  of  simulations,  all  of  which  were  sample  size  n  =  1. 

Deviations  from  the  normal  approximation  for  some  parameters  was  indicated  by  the 
Monte  Carlo  simulation  results.  As  was  observed  from  the  histograms,  the  concentrations 
of  Na  and  alkalinity  in  Chub  Lake  were  negatively  skewed  as  all  or  most  of  the  counts  were 
above  the  upper  confidence  limit  (%>  +L).  CI  concentrations  during  ice  cover  were  also 
negatively  skewed  and  pH  was  positively  skewed  in  Blue  Chalk  Lake.  An  apparent 
difference  between  seasons  was  observed  in  Harp  Lake  for  colour.  The  counts  suggest 
colour  was  positively  skewed  during  the  SS  period  and  negatively  skewed  during  the  IC 
period.  However,  in  all  the  above  cases  the  percentage  of  counts  remained  below  or  near 
5%. 

Large  deviations  from  the  expected  normal  distribution  are  indicated  for  DOC 
concentrations  in  Harp  Lake  and  Al  concentrations  in  all  three  lakes.  DOC  concentrations 
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were  negatively  skewed,  as  most  counts  were  beyond  the  upper  95%  confidence  interval, 
with  an  overall  estimated  a  of  0.056  to  0.083.  Similarly,  Al  concentrations  were  negatively 
skewed  in  all  three  lakes,  very  strongly  in  Blue  Chalk,  with  a  relatively  large  %  of  the 
counts  outside  the  95%  C.I.  (4.3  to  8.5%). 

iii)       Sample  Size 

1)         The  Normal  Approximation 

The  numbers  of  samples  required  to  estimate  a  submean  which  is  within  a  given 
distance  of  the  long  term  mean,  with  95%  confidence,  for  the  3  study  lakes  (calculated  from 
1976  to  1984)  are  shown  in  Tables  7  and  8.  The  95%  C.I.  are  expressed  as  a  function  of 
sample  size  for  the  two  seasons:  summer  stratified  (SS)  and  ice  cover  (IC)  in  Fig.  7. 

For  most  parameters,  one  random  sample  would  probably  be  within  20  or  30%  of 
the  annual  mean  for  the  study  lakes.  Similarly,  sampling  frequencies  of  1  to  3  per  year 
would  produce  mean  annual  estimates  within  20%  of  the  long  term  mean  about  95%  of  the 
time.  The  two  exceptions  are  the  parameters  colour  and  Al  which  have  large  variability 
about  the  mean  and,  thus,  require  considerably  greater  sampling  effort.  To  be  within  10% 
of  the  long  term  mean,  an  average  calculated  from  up  to  20  samples  may  be  required  for 
the  12  parameters.  If  a  larger  confidence  interval  than  95%  is  acceptable,  then  the  number 
of  samples  can  be  reduced.  The  variability  associated  with  seasonal  means  are  comparable 
to  the  long  term  annual  means  and  the  similarity  in  estimated  sampling  frequenc>'  for  the 
given  confidence  interval  (95%)  reflects  this.  Analytical  and  tiualitv  ciMitrol  standard 
deviations  are  so  small  relative  to  (lie  natural  variability  that  intensive  sainfilink:  treiiueneie.s 
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would  be  required  to  produce  estimates  as  close  as  the  analytical  and  quality  control 
standard  deviations  from  the  long-term  mean. 

The  standard  deviations  for  conductivity  are  very  consistent  (range  0.76  to  2.00)  for 
the  three  lakes  and  two  seasons.  One  sample,  taken  any  time  of  the  year,  would  be  within 
2  to  4  nS,  or  about  10%,  of  the  long  term  annual  mean. 

Similarly  the  standard  deviation  for  alkalinity  are  consistently  around  0.2  mg  1'^  and 
1  random  sample  would  be  within  0.4  to  0.75  mg  1'^  of  the  annual  mean.  A  sampling 
frequency  of  2-4  per  year  would  produce  estimates  within  10%  of  the  long  term  aimual 
mean.  However,  the  annual  mean  for  Chub  Lake  is  considerably  lower.  Thus,  for  dilute 
Chub  Lake  a  sampling  frequency  of  approximately  10  would  be  required  to  be  within  20% 
of  the  annual  mean. 

The  annual  and  seasonal  variability  of  pH  is  consistent  between  lakes  with  a  standard 
deviation  ranging  from  0.139  to  0.198  pH  units.  One  random  sample  would  be  within  0.4 
units,  or  about  5%,  of  the  aimual  mean.  A  more  acceptable  tolerance  of  0.1  units  from  the 
mean  would  require  as  many  as  8  to  16  samples  per  year  or  season.  A  sampling  frequency 
of  2  to  4  should  estimate  an  aimual  or  seasonal  mean  within  0.2  imits  of  the  actual  mean, 
95%  of  the  time. 

The  standard  deviation  for  seasonal  and  annual  means  for  Ca  concentration  ranged 
from  0.247  to  0.389.  To  get  an  estimate  of  seasonal  or  annual  Ca  concentration  within  10% 
(0.2  to  0.4  mg  r^)  of  the  mean,  95%  of  the  time  you  would  require  2  to  7  samples.   One 
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or  2  random  samples  would  produce  an  estimate  within  0.5  to  0.7  mg  r\  or  about  20%,  of 
the  mean.   Chub  has  a  more  variable  Ca  distribution  in  both  IC  and  SS  periods. 

The  variation  for  Mg  is  very  consistent  (standard  deviation  0.072  to  0.085  mg  1"^) 
and  a  sampling  frequency  of  2  to  5  would  produce  a  seasonal  or  annual  estimate  within 
10%  of  the  actual  mean.  One  random  sample  would  produce  an  estimate  within  about 
0.15  mg  r^  (20  to  25%)  of  the  actual  seasonal  or  annual  mean.  Similarly,  potassium  has 
a  consistent  variability  for  the  lakes  and  seasons  considered  here  (0.044  to  0.066  mg  K  f 
^).  Six  samples  are  likely  to  estimate  within  the  analytical  standard  deviation  (0.04mg  1'^) 
of  the  seasonal  or  annual  mean.  One  or  2  samples  would  estimate  within  0.1  mg  V'  (20  to 
30%)  of  the  annual  mean. 

The  standard  deviation  of  the  annual  means  for  Na  was  consistent  between  lakes 
but  there  appears  to  be  slightly  more  variability  during  IC  than  SS  seasons,  except  in  Blue 
Chalk  Lake  where  it  is  the  same  year  round.  For  annual  means,  approximately  4  to  13 
samples  and  2  to  3  samples  would  be  required  to  estimate  within  10%  (0.07  to  0.10  mg 
r^)  and  20%  of  the  actual  mean,  respectively.  One  random  sample  would  be  within  about 
0.25  mg  r^  (20  to  35%)  of  the  actual  annual  mean.  Chub  Lake  during  IC  is  the  most 
variable  example  and  a  slightly  greater  sampling  frequency  of  about  17  and  4  would  be 
required  to  estimate  within  10%  (0.08  mg  V^)  and  20%  (0.15  mg  1"^),  respectively,  of  the 
actual  mean. 

Chloride  standard  deviations  are  relatively  large  compared  to  the  means  and  ranged 
from  0.060  to  0.131  nig  1  '.    A  sampling  frec|uency  ot  ^'  to  25  would  be  needed  to  eslinuiie 
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within  10%  of  the  actual  annual  mean.  Two  to  6  samples  would  estimate  within  20%  (0.09 
to  0.25  mg  r^)  of  the  actual  mean  95%  of  the  time.  As  with  Na,  there  appears  to  be  more 
variability  during  IC  than  SS,  except  in  Harp  Lake  where  it  is  the  same  year  round.  Chub 
Lake  during  IC  is  the  most  variable  example. 

The  standard  deviation  for  SO4  was  relatively  consistent  between  lakes  and  seasons, 
ranging  from  0.292  to  0.474  mg  1"^  for  data  analyzed  by  ion  chromatography  only.  The 
variance  is  low  relative  to  the  long-term  mean  and  1  or  2  samples  would  estimate  within 
10%  of  the  actual  aimual  or  seasonal  mean. 

There  is  considerable  difference  in  DOC  variability  between  the  three  lakes.  Blue 
Chalk  and  Chub  Lakes  had  lower  standard  deviation  and  approximately  3  to  7  samples 
would  estimate  within  10%  of  the  actual  annual  mean.  Just  1  or  2  samples  would  be 
required  to  estimate  within  20%  of  the  annual  mean.  DOC  concentrations  in  Harp  Lake 
were  much  more  variable  relative  to  the  mean.  This  was  primarily  due  to  high  variability 
during  SS  (0.887  mg  1'^)  compared  to  IC  (0.341  mg  V^).  To  estimate  within  10%  of  the 
actual  seasonal  mean  DOC  in  Harp  Lake,  18  samples  would  be  required  during  SS 
compared  to  2  or  3  during  IC.  Just  1  sample  during  IC  would  be  within  20%  of  the  actual 
seasonal  mean,  while  it  would  take  3  or  4  samples  during  SS  to  get  similar  precision. 

Both  colour  and  Al  have  a  large  standard  deviation  relative  to  the  annual  and 
seasonal  mean  in  all  three  study  lakes.  An  intensive  sampling  regime  would  be  required 
to  have  acceptable  precision  on  estimates  of  average  concentrations.  Greater  than  10 
samples  for  colour  and  30  samples  for  Al  would  be  required  to  estimate  within  10%  of  the 
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annual  or  seasonal  mean.  One  random  sample  would  be  greater  than  35%  and  50%  of  the 
actual  mean  for  colour  and  Al,  respectively.  The  relative  variability  does  decrease  as  the 
colour  or  Al  concentration  increases.  In  highly  colored  Chub  Lake  approximately  3  to  5 
samples  would  estimate  within  20%  of  the  actual  mean.  Chub  Lake  also  had  the  highest 
Al  concentration,  but  at  least  7  samples  would  be  required  to  estimate  within  20%  of  the 
actual  annual  or  seasonal  mean. 

2)        Non-parametric  Estimates 

An  attempt  was  made  to  estimate  the  lowest  number  of  samples  required  based  on 
the  95%  quantile  and  the  binomial  distribution  (Cotter  1985).  This  method  has  the 
advantage  of  being  distribution  free.  However,  the  method  requires  a  minimum  of  5  to  6 
samples  to  define  a  confidence  limit  about  a  quantile  and  that  appears  to  be  a  larger 
number  of  samples  than  required  as  indicated  by  the  data. 


DISCUSSION 

i)  Assumptions 

The  practical  application  of  these  results  require  acknowledgement  of  both  ihe 
statistical  assumptions  inherent  in  the  methodology  and  of  what  we  can  call  the  assumptions 
resulting  from  regionalization  or  extrapolation  of  the  findings. 
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The  statistical  methods  employed  here  assume  1)  the  data  are  independent  and 
random  samples  of  the  true  population,  and  2)  for  some  of  the  statistical  tests,  that  the 
measured  values  for  each  parameter  are  normally  distributed;  thus,  the  parametric  statistics 
can  be  used  to  estimate  u,  variance  and  the  sample  size  required. 

As  noted  earher,  the  lakes  were  not  randomly  sampled.  In  practical  terms  a 
completely  random  data  set  is  extremely  difficult,  if  not  impossible,  to  obtain. 
Systematically  sampled  data  can  be  treated  as  random  sometimes  as  long  as  there  is  no 
systematic  behaviour  related  to  the  parameters  measured  (Ellis  and  Lacey  1980).  If  there 
is  a  systematic  trend,  then  it  is  important  that  the  sampling  regime  be  out  of  step  with  any 
naturally  occurring  frequency  in  the  lakes  if  the  risk  of  significant  bias  is  to  be  avoided. 
However,  it  is  imlikely  that  natural  variation  would  be  in  cycle  with  a  regular  'workweek'. 
In  addition  the  sequence  of  sampling  varied  year  to  year  and  the  irregularity  of  natural 
factors  would  reduce  if  not  prevent  alignment  with  seasonal  or  perhaps  shorter  cycles.  Thus 
the  multi-year  pooled  data  sets  are  assumed  to  represent  an  intensive  random  sampling  of 
individual  lakes  for  an  average  hydrologic  year. 

The  danger  or  bias  of  systematic  sampling  when  sampling  frequency  and  natural 
cycling  coincide  is  that  it  may  artificially  reduce  estimates  of  true  variance  in  the  population. 
However,  in  this  study  artificially  reduced  variance  would  result  in  a  predicted  lower 
sampling  frequency  needed  for  a  given  95%  confidence  interval  about  the  mean  assuming 
random  sampling.  Just  how  liberal  an  estimate  of  required  sample  size  depends  on  the  size 
of  the  cycle  in  relation  to  the  variability  of  the  time  interval   (i.e.  1  or  2  per  month).   As 
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mentioned  above,  we  have  assumed  that  the  sampling  frequency  does  not  coincide  with 
naturally  occurring  cycles. 

Serial  dependence  of  the  data  is  a  common  and  can  be  a  serious  problem  in  water 
quality  surveys.  Generally,  temporal  correlations  in  water  quality  measures  are  sHght  if  the 
time  interval  is  relatively  large  (Millard  et  al.  1985).  However,  SS  sampling  frequencies  in 
this  study  were  often  greater  than  1  per  month  and  potential  serial  dependence  may  prevent 
the  application  of  a  distribution  model  (Lettenmaier  1976). 

If  we  believe  that  the  multi-year  pooled  data  set  represents  random  samples  about 
a  relatively  stable  long-term  mean,  then  we  have  a  data  set  large  enough  to  adequately 
determine  sample  distribution.  The  results  from  the  Monte  Carlo  simulation  suggest  that 
most  parameters  approximate  a  normal  distribution.  It  appears  that  the  parametric 
methods  employed  here  are  robust  enough  to  accommodate  some  departure  from  normality. 
Other  workers  have  suggested  that  assuming  normal  distribution  is  valid  in  water  qualit>' 
data  as  long  as  one  is  concerned  with  arithmetic  mean  (Montgomery  and  Hart  1974,  Walker 
1982).  This  is,  in  part,  a  statement  of  the  central  limit  theorem  that  whatever  the  shape  of 
the  frequency  distribution  of  the  original  population,  the  distribution  of  means  in  repeated 
random  samples  of  size  'n'  tends  to  become  normal  as  'n'  increases.  It  would  seem  that 
even  a  low  'n'  of  2  or  3  will  suffice  to  make  the  2.5  and  97.5  percentiles  of  the  lake 
chemistry  means  fall  on  or  near  the  normal  95%  confidence  interval  boundaries. 

Artificial  truncation  of  the  distribution  by  a  detection  linnt  may  introduce 
considerable   bias.     This   may   be   a   problem  with   sonic   ot    the   exaniiiiod   parameters. 
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especially  Al  (i.e.  Blue  Chalk  Lake  data).  Monte  Carlo  simulations  suggest  that  Al  and 
possibly  DOC  data  deviate  substantially  from  the  normal  distribution  and  that  perhaps  a 
log-normal  transformation  should  be  performed  or  non-parametric  statistics  used  (Cotter 
1985).  Care  should  be  taken  with  any  case  where  a  mean  concentration  is  near  the 
detection  limit  of  the  analytic  method. 

Regionalization  or  extrapolation  assumptions  concerning  the  sampling  procedure 
include: 

1)  the  best  estimate  of  the  true  value  for  the  concentration  of  the  12  parameters 
is  our  volume-  weighted  annual  average  for  our  intensive  study  lakes, 

2)  individual  samples  are  a  true  reflection  of  water  quality  at  time  of  sampling, 

3)  there  is  no  difference  in  true  values  of  the  12  parameters  between  years  for 
Blue  Chalk,  Harp  and  Chub  Lakes  and,  therefore,  the  multiple  year  data  set 
can  be  lumped  to  increase  sample  size, 

4)  the  three  intensively  studied  lakes  represent  a  large  number  of  lakes  in  this 
province  and  results  from  them  may  be  extrapolated  to  some  range  of  lakes 
with  similar  chemistry  and  can  be  useful  in  determining  survey  sampling 
strategies, 

5)  The  discussion  only  applies  to  our  particular  field  methods.  Other  field 
methods  may  not  have  the  same  variance  je.  grab  samples  vs.  volume 
weighted  samples. 
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The  representativeness  of  a  point  sample  is  an  important  consideration  in  any  study 
but  is  beyond  the  scope  of  this  manuscript.  Ellis  and  Lacey  (1980),  Walker  (1982)  and  Mar 
.et  al  (1986)  give  relevant  discussions  concerning  this  topic. 

Lumping  the  multi-year  data  set  was  done  to  increase  the  sample  size  in  order  to 
more  effectively  determine  the  distribution  of  the  data  and  population  variance.  Ideally, 
approximately  100  samples  should  have  been  taken  randomly  from  each  lake  during  one 
hydrologic  year.  With  multiple  years  it  may  be  difficult  to  distinguish  between  trend  and 
random  variation.  The  three  lakes  were  deliberately  chosen  because  there  were  little  or 
no  annual  trends.  However,  the  difference  between  years  will  simply  increase  our  estimate 
of  variability  and  any  conclusion  drawn  from  the  estimate  will  tend  to  be  conservative  when 
applied  as  a  within  year  variability. 

The  representative  nature  of  the  3  studied  lakes  to  the  population  of  lakes  to  which 
the  results  can  be  applied  is  a  very  important  consideration  if  the  results  are  to  be 
extrapolated  to  other  less  intensively  studied  water  bodies.  The  mean  annual  concentration 
of  the  12  parameters  in  the  three  lakes  roughly  cover  the  range  found  in  most  lakes  in  the 
area.  Similar  physiographic  (including  climatic  variability)  conditions  occur  over  the  area 
and  presumably  most  lakes  will  have  similar  within  year  and  between-year  variability.  Of 
the  lakes  in  the  study  area  that  were  sampled  frequently  enough  to  estimate  variance,  most 
have  variances  within  the  range  of  the  three  study  lakes.  These  results  may  also  be 
extrapolated,  with  some  discretion,  to  other  water  bodies  with  similar  mean  annual 
concentrations,  geologic  and  physiographic  characteristics. 
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ii)         General  Remarks 

To  determine  the  sampling  frequency  for  a  given  sampling  programme  on  a  lake  the 
following  information  is  required: 

1)  how  much  uncertainty  can  be  tolerated,  and  to  what  degree  of  confidence,  in 
the  estimate,  and 

2)  what  is  the  variabiUty  of  the  water  chemistry. 

It  appears  that  most  of  the  examined  parameters  are  roughly  normally  distributed 
about  the  mean  annual  concentration.  If  enough  data  are  present  to  calculate  past 
estimates  of  variability  of  a  lake,  then  parametric  methods,  such  as  employed  here,  may  be 
employed  to  determine  the  sampling  frequency  required  with  reasonable  safety.  Problems 
of  systematic  sampling  and  serial  dependence  of  the  data  must  still  be  considered.  The 
exceptions  are  Al  and  possibly  DOC,  which  appear  to  have  somewhat  skewed  distributions 
in  the  three  lakes  studied.  However,  the  variance  was  considerably  high  for  these  two 
parameters  and  estimates  of  sampling  frequency  were  high  as  well.  Considerable  skewness 
may  result  from  truncating  when  lake  water  concentrations  are  near  detection  limits.  In 
these  situations  a  log-normal  approximation  or  non-parametric  methods  should  be 
employed. 

If  adequate  information  does  not  exist  about  the  variance  of  a  given  lake  chemistry 
then  the  graphs  may  be  used  to  estimate  the  sampling  frequency  required  for  a  given  survey 
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program  to  produce  estimates  within  a  given  tolerance  of  the  true  mean  with  95% 
confidence.  The  graphs  depicting  the  95%  confidence  interval  as  a  function  of  sample  size 
show  that  the  estimates  rapidly  converge  to  the  true  mean  from  sample  sizes  one  to  four, 
with  little  improvement  beyond  that.  Provided  the  lake  has  similar  chemistry  and 
physiographic  characteristics  to  the  study  lakes  these  should  be  reasonable  estimates. 

Marked  seasonal  trends  occur  in  most  lake  water  parameters  measured  in  lakes. 
The  results  presented  here  suggest  that  there  is  very  little  difference  between  various 
seasons  and  annual  means  on  the  long  term.  Ellis  and  Lacey  (1980)  suggest  expressing  the 
systematic  time-dependence  of  concentration  as  a  ratio  (within  period  standard  deviation) 
/  (overall  standard  deviation).  The  smaller  the  number  the  greater  the  '  benefit'  or  bias 
of  systematic  sampling.  The  standard  deviation  of  the  two  seasons  with  ample  data,  SS  and 
IC,  are  as  large  or  larger  than  the  annual  variability  suggesting  little  bias  on  variance 
estimates  in  systematically  sampling  at  least  on  a  seasonal  basis  for  these  lakes.  This  is  true 
of  conductivity,  alkalinity,  Ca,  Mg,  K,  and  SO,  where  1  or  2  samples  would  be  within  20%, 
and  often  10%,  of  the  aimual  mean.  Some  parameters  have  consistently  lower  or  higher 
concentrations,  but  not  significantly  so,  during  certain  seasons.  This  information  would  be 
useful  in  determining  when  to  sample  for  maximum  or  minimum  annual  values. 

The  annual  variability  for  DOC,  Colour  and  A!  was  considerable  and  large  sampling 
frequencies  are  required  to  estimate  annual  and  seasonal  means  with  reasonable  precision. 
This  is  al.so  the  case  for  pH  if  an  acceptable  tolerance  of  0.1  or  0.2  pH  units  is  considered. 
Systematic  sampling  may  be  employed  to  separate  the  une.vplained  variation  and  expliMi 
the  regularly  occurring  seasonal  cycles  in  water  chonlisl^^■  in  the  sampling  design  (I'llis  and 


27 


Lacey  1980).  It  appears  that  the  period  of  the  within-year  or  seasonal  cycles  is  shorter  then 
the  two  defined  hydrologic  seasons  (SS  and  IC)  and  each  season  incorporates  a  large 
portion  of  the  overall  annual  variance.  Utilizing  more  detailed  information  of  the  within- 
year  cycles  in  a  systematic  sampling  regime  could  result  in  substantial  reduction  in  sampling 
error  and  sample  size  needed  for  a  given  precision. 

Seasonal  trends  in  water  chemistry  may  be  more  pronounced  in  smaller  or  more 
sensitive  waterbodies,  possibly  due  to  the  greater  interaction  with  seasonal  cycles  of  runoff 
chemistry  with  the  surrounding  watershed.  Chub  Lake,  which  is  the  most  dilute  and 
potentially  the  most  sensitive  to  acid  deposition,  consistently  had  larger  variance  relative 
to  the  annual  mean.  Seasonal  means  were  also  consistently  different  from  the  annual 
means,  especially  during  IC  periods.  Greater  sampling  frequency  of  such  lakes  or 
exploitation  of  seasonal  cycles  to  reduce  the  required  sample  size  may  be  necessary  for 
surveys  of  potentially  sensitive  lakes. 
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Table  1.       Location  and  morphometric  features  of  the  study  lakes. 


Blue  Chalk  Harp  Chub 


Latitude/Longitude 

45' 

'l2'N  78°56'W 

45°23'N  79°08'W 

45" 

'l3'N78°59'W 

Lake  Area  (km^) 

0.494 

0.669 

0.322 

Watershed  Area  (km^) 

1.35 

5.09 

2.86 

Mean  Depth  (m) 

8.5 

12.4 

8.9 

Maximum  Depth  (m) 

23 

40 

27 

Table  2.       Chemical  parameters,  standard  deviation  of  the  test  and  analytical  method  used  (Ontario 
Ministry  of  the  Environment  1983,  1985). 


Parameter 

Standard  Deviation 

Method 

Cond 

0.22   ^iS 

radiometry  meter 

Alkalinity 

0.050  mg  L-' 

CaC03 

inflection  point  potentiometer  titration 

PH 

0.032  pH 

electrometric 

Ca 

0.151  mg  L-' 

atomic  absorption 

Mg 

0.042  mg  L-= 

atomic  absorption 

Na 

0.070  mg  L-' 

atomic  absorption 

K 

0.040  mg  L-' 

atomic  absorption 

CI 

0.051  mgL-' 

colourimetry 

SO, 

0.536  mg  L-' 

1976-1980  methylthymol  blue  colourimetry 

0.167  mg  L-' 

1980-1985  ion  chromatography 

DOC 

0.148  mgL"' 

UV  digestion,  colourimetry 

Colour 

0.52   Hazen  units 

colourimetry 

Al  (total) 

5.3     /igL-' 

atomic  adsorption 

Table  3.  Mean  values  and  standard  deviations  for  12  chemical  parameters  in  three  lakes  in  the 
annual  (AN),  summer  stratified  (SS)  and  ice-covered  (IC)  seasons.  Values  include  data 
collected  from  1976  to  1985. 


Lake 

AN 

SS 

IC 

Parameter 

N 

mean 

s.d. 

N 

mean 

s.d. 

N 

mean 

s.d. 

Cond 

BC 

238 

29.9 

1.35 

179 

29.9 

1.40 

38 

30.4 

1.16 

HP 

247 

35.3 

1.11 

176 

35.2 

0.76 

42 

35.9 

1.81 

CB 

244 

29.7 

1.82 

178 

29.6 

1.74 

42 

30.6 

2.00 

Alkalinity 

BC 

120 

3.99 

0.369 

87 

4.00 

0.39 

19 

4.06 

0.18 

HP 

122 

3.21 

0.225 

87 

3.19 

0.21 

19 

3.28 

0.22 

CB 

116 

0.70 

0.212 

84 

0.70 

0.22 

20 

0.70 

0.21 

PH 

BC 

238 

6.59 

0.201 

179 

6.64 

0.198 

38 

6.42 

0.147 

HP 

248 

6.30 

0.145 

177 

6.32 

0.142 

42 

6.25 

0.139 

CB 

244 

5.66 

0.179 

179 

5.68 

0.182 

42 

5.62 

0.183 

Ca 

BC 

201 

2.76 

0.274 

147 

2.77 

0.286 

34 

2.74 

0.252 

HP 

196 

3.12 

0.260 

135 

3.11 

0.247 

33 

3.18 

0.309 

CB 

199 

2.68 

0.359 

138 

2.67 

0.356 

38 

2.76 

0.389 

Mg 

BC 

200 

0.74 

0.079 

146 

0.74 

0.081 

34 

0.75 

0.067 

HP 

200 

0.93 

0.080 

135 

0.94 

0.078 

37 

0.95 

0.085 

CB 

197 

0.68 

0.074 

135 

0.68 

0.078 

39 

0.70 

0.072 

Na 

BC 

201 

0.78 

0.210 

148 

0.77 

0.12 

33 

0.83 

0.11 

HP 

202 

1.00 

0.111 

135 

0.98 

0.10 

39 

1.06 

0.12 

CB 

202 

0.70 

0.124 

138 

0.68 

0.11 

41 

0.76 

0.16 

K 

BC 

201 

0.37 

0.05 

149 

0.37 

0.05 

32 

0.38 

0.04 

HP 

201 

0.54 

0.05 

135 

0.53 

0.04 

39 

0.55 

0.07 

CB 

202 

0.40 

0.05 

138 

0.40 

0.05 

41 

0.41 

0.06 

CI 

BC 

204 

0.43 

0.065 

149 

0.42 

0.060 

36 

0.46 

0.077 

HP 

203 

0.77 

0.118 

135 

0.77 

0.120 

40 

0.78 

0.124 

CB 

201 

0.49 

0.122 

137 

0.49 

0,117 

42 

0.51 

0,131 

SO, 

BC 

94 

6.79 

0.453 

69 

6.79 

0.474 

15 

6.93 

0.374 

HP 

93 

8.08 

0.380 

65 

8.10 

0.410 

15 

8.10 

0.292 

CB 

93 

7.56 

0.396 

64 

7.52 

0.393 

18 

7.74 

0.423 

DOC 

BC 

137 

1.8 

0.227 

99 

1.9 

0.208 

23 

1.8 

0.182 

HP 

144 

4.1 

0.856 

96 

4.1 

0.887 

26 

4,1 

0,341 

CB 

129 

4.9 

0.473 

88 

4.8 

0,430 

25 

5  1 

0450 

Colour 

BC 

91 

6.0 

2.37 

64 

6.8 

\MA 

17 

3  6 

1  87 

HP 

95 

21.1 

4.46 

65 

20.2 

4.21 

17 

23  1 

3  4.S 

CB 

90 

46.5 

8.63 

62 

45.1 

8.01 

18 

4^,5 

»)  16 

OAI 

HC 

^M 

1  1 .62 

10.80 

63 

10, .S3 

8,33 

1^ 

1  \  ^> 

1 :  ^8 

MP 

107 

43.l'> 

15  43 

68 

4:(. 

\^^  :-» 

:3 

400 

14  "4 

CB 

111 

95.'>: 

20,60 

73 

'>.V8 

:-i,:s 

:.i 

i)>.i) 
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Table  4.  Comparison  of  seasonal  mean  values  to  the  hydrologic  year  mean  and  test  of  significant 
differences.  The  nonparametric  test  used  is  the  Wilcoxon  matched-pairs  signed-ranks 
test.  Results  are  reported  as  the  means  of  the  differences  (seasonal  means-annual  mean) 
and  the  significance  level  of  the  Wilcoxon  test.  (SO  =  spring  overturn;  SS  =  summer 
stratified;  FO  =  fall  overturn;  IC  =  ice-covered;  IF  =  ice-free;  AN  =  annual). 


Parameter 

Lake 

SO-AN 

SS-AN 

FO-AN 

IC-AN 

IF-AN 

Cond 

BC 
HP 
CB 

-0.7* 

-0.6 

-0.5 

-0.1 
-0.0 
-0.1 

-0.4 

-0.3 

0.4 

0.4 
0.3 
0.7** 

-0.1 

0.0 

-0.1 

Alkalinity 

BC 
HP 
CB 

-2.81 
-0.19 
-0.12 

-0.11 

0.01 

-0.02 

0.06 
0.21* 

0.22* 

0.14 

0.07 

-0.00 

0.02 
-0.01 

pH 

BC 
HP 
CB 

-0.06 
-  0.06* 
-0.04 

0.04* 

0.02 

0.02 

-0.05 
-0.02 
0.09* 

-0.17* 

-0.04 

-0.04 

0.02 
0.01 
0.02 

Ca 

BC 
HP 
CB 

-0.17* 

-0.04 

-0.04 

0.01 
-0.02 
-0.01 

0.11 

0.06 

-0.06 

-0.04 
0.04 
0.08 

0.01 

-0.01 

0.02 

Mg 

BC 
HP 
CB 

-  0.05* 

-  0.05* 
0.02 

0.00 

0.00 

-0.00 

-0.04 

-0.01 

0.02 

0.01 

0.02* 

0.02** 

-0.00 

-0.00 

0.00 

Na 

BC 
HP 
CB 

0.03 
0.02 
0.01 

-  0.01 

-  0.02* 

-  0.02* 

0.06 

-0.03 

0.09 

0.04** 

0.06 

0.05* 

-0.01 

-0.02 

0.01 

K 

BC 
HP 
CB 

-  0.003 

-  0.003 

-  0.002 

0.002 

-  0.003 

-  0.001 

0.01 

-0.01 

0.02 

0.01 
0.01 
0.01 

-0.0003 

-  0.002 

0.001 

CI 

BC 
HP 
CB 

-0.03 

0.01 

-0.03 

0.004 
-0.01 
-  0.0026 

O.OI 

0.01 

-  0.02 

0.03* 

0.01 

0.01 

-  0.004 
-0.01* 

-  0.001 

so« 

BC 
HP 
CB 

-0.10 
-0.13 
-0.15 

-0.01 

0.001 
-0.07 

0.10 

-0.02 

0.02 

0.10 
0.17 
0.27* 

-0.02 
-  0.002 
-0.08 

DOC 

BC 
HP 
CB 

-0.12 
-0.21 
-  0.33* 

-0.03 
-0.03 
-0.05 

-0.05 
0.09 
0.16 

-0.08 
-0.03 
0.24* 

-0.05 
-0.05 
-0.05 

Colour 

BC 
HP 
CB 

-0.4 

0.2 

-2.6 

0.3 

-0.8 

2.3 

1.0 
2.4 
4.4 

-2.3 
2.0* 
0.7 

-0.3 

-0.5 

2.7 

Al 

BC 
HP 
CB 

2.2 
1.1 
7.6 

-0.08 
0.4 
-0.3 

2.9 
-0.9 

-  2.2 

3.6 
5.6 

3.1 

2.2 
0.3 

3.9 

* 
** 

p  <  0.05 
p  <  0.01 

Table  5.  Levels  of  significance  of  the  results  from  the  Lilliefors  test  of  normality.  Not  significant 
(n.s.)  indicates  the  data  fit  the  normal  distribution,  and  one  or  two  stars  indicate  the 
Lilliefor's  statistic  exceeds  the  0.95  or  0.99  quantile  respectively. 


Parameter 


Blue  Chalk 


Harp 


Chub 


Annual 


Conductivity 

** 

Alkalinity 

* 

PH 

** 

Ca 

* 

Mg 

* 

Na 

** 

K 

** 

CI 

** 

SO, 

n.s 

DOC 

n.s, 

Colour 

n.s, 

Al 

** 

Summer  Stratified 

Conductivity 

** 

Alkalinity 

n.s. 

pH 

*• 

Ca 

« 

Mg 

•* 

Na 

* 

K 

** 

CI 

n.s. 

SO, 

•« 

DOC 

•* 

Colour 

** 

Al 

** 

Ic?  C9v?r 

Conductivity 

Alkalinity 

pH 

Ca 

Mg 

Na 

K 

CI 

SO, 

DOC 

Colour 

Al 


n.s. 
n.s. 
n.s. 
n.s. 
n.s. 
n.s. 

n.s. 

•• 

n.s. 
n.s. 
n.s. 


n.s. 

** 

** 
** 
** 
** 

** 
*« 

n.s. 


n.s. 


** 
** 
** 

n.s. 

«* 

n.s. 


n.s. 
n.s. 
n.s. 
n.s. 
n.s. 
n.s. 
n.s. 
n.s. 
n.s. 
n.s. 


** 
** 
** 
*« 
** 
«* 
«* 

n.s. 

n.s. 

«• 


n.s. 

* 

n.s. 
n.s. 


n.s. 


n.s. 
n.s. 
n.s. 
n.s. 
n.s. 
n.s. 
n.s. 

n.s. 
n.s. 
n.s. 
n.s. 


Table  6a.  Comparison  of  Monte  Carlo  simulation  with  the  parametric  confidence  limit  for 
conductivity.  L  is  the  97.5%  interval  about  the  mean  assuming  a  normal  distribution. 
Number  of  cases  in  the  simulation  above  and  below  beyond  the  confidence  interval  were 
counted  and  expressed  as  a  %.  a  is  the  probability  of  finding  sample  mean  outside  the 
95%  confidence  interval. 


Lake 


Season' 


Sample 

Number  of 

Monte  Carlo 

Size 

Cases 

L 

%<(x-L) 

%>(x+L) 

Q 

1 

247 

7.255 

0 

0.4 

0.004 

2 

30,381 

5.130 

0 

0.8 

0.008 

3 

2,481,115 

4.189 

0 

1.2 

0.012 

1 

179 

2.741 

0.6 

3.4 

0.039 

2 

15,931 

1.938 

1.1 

3.3 

0.045 

3 

939,929 

1.583 

1.7 

3.5 

0.051 

1 

38 

2.283 

2.6 

2.5 

0.053 

2 

703 

1.614 

3.7 

5.0 

0.087 

3 

8,436 

1.318 

2.0 

2.3 

0.043 

1 

244 

3.559 

0.8 

1.6 

0.025 

2 

29,646 

2.517 

1.9 

2.3 

0.042 

3 

2,391,444 

2.055 

2.4 

2.4 

0.048 

1 

176 

8.352 

0 

0.6 

0.006 

2 

15,400 

5.906 

0 

1.1 

0.011 

3 

893,200 

4.822 

0 

1.7 

0.017 

1 

42 

3.550 

0 

4.8 

0.048 

2 

861 

2.510 

0.5 

6.6 

0.071 

3 

11,480 

2.050 

0.7 

5.4 

0.061 

1 

238 

2.645 

0.8 

2.9 

0.038 

2 

28,203 

1.870 

1.1 

3.4 

0.045 

3 

2,218,636 

1.527 

1.6 

3.3 

0.049 

1 

178 

3.414 

1.1 

1.1 

0.022 

2 

15,753 

2.414 

2.3 

1.8 

0.041 

3 

924,176 

1.971 

2.6 

2.0 

0.046 

1 

42 

3.926 

0 

2.4 

0.024 

2 

861 

2.776 

1.5 

3.2 

0.047 

3 

11,480 

2.267 

1.6 

3.0 

0.046 

Blue  Chalk 


Harp 


Chub 


AN 


SS 


IC 


AN 


SS 


IC 


AN 


SS 


IC 


AN  =  annual;  SS  =  summer  stratification;  IC  =  ice-covered. 


Table  6b.     Comparison  of  Monte  Carlo  simulation  count  of  number  of  Monte  Carlo  cases  beyond 
the  95%  confidence  limit  (L)  for  alkalinity. 


Sample 

Number  of 

Monte  Carlo 

Lake 

Season' 

Size 

Cases 

L 

%<(x-L) 

%>(x+L) 

Or 

Blue  Chalk 

AN 

1 

120 

0.748 

1.7 

1.7 

0.033 

2 

7,140 

0.529 

2.4 

2.0 

0.044 

3 

280,840 

0.432 

2.5 

2.2 

0.047 

SS 

1 

87 

0.799 

2.3 

2.3 

0.046 

2 

3,741 

0.565 

2.3 

2.0 

0.043 

3 

105,995 

0.461 

2.4 

2.2 

0.046 

IC 

1 

19 

0.358 

0 

0 

0 

2 

171 

0.253 

2.8 

0 

0.028 

3 

969 

0.207 

2.9 

0.4 

0.032 

Harp 

AN 

1 

122 

0.441 

0.8 

2.5 

0.033 

2 

7,381 

0.312 

0.9 

3.4 

0.043 

3 

295,240 

0.254 

1.2 

3.6 

0.048 

SS 

1 

87 

0.413 

1.1 

1.1 

0.023 

2 

3,741 

0.292 

1.0 

2.6 

0.036 

3 

105,995 

0.238 

1.2 

3.5 

0.047 

IC 

1 

19 

0.438 

0 

5.3 

0.053 

2 

171 

0.310 

2.2 

4.2 

0.064 

3 

969 

0.253 

1.3 

3.1 

0.044 

Chub 

AN 

1 

116 

0.416 

0 

4.3 

0.043 

2 

6,670 

0.294 

0.6 

4.3 

0.048 

3 

253.460 

0.240 

0.8 

4.3 

0.051 

SS 

1 

84 

0.421 

0 

3.6 

0.036 

2 

3,486 

0.298 

0.6 

4.5 

0.052 

3 

95.284 

0.243 

0.7 

4.8 

0.055 

IC 

I 

20 

0.409 

0 

5.0 

0.050 

2 

190 

0.289 

0 

4.0 

0040 

3 

1.140 

0.23') 

1,0 

3.1 

0  041 

AN  -  annual;  SS  -  summer  stratification;  IC  -  ice-covered. 


Table  6c.     Comparison  of  Monte  Carlo  simulation  with  the  parametric  confidence  limit.    Number 
of  cases  in  the  simulation  beyond  the  95%  confidence  interval  were  counted  for  pH. 


Lake 


Season' 


Sample 
Size 


Number  of 
Cases 


Monte  Carlo 
%<(x-L)      %>(x+L) 


Blue  Chalk 


Harp 


Chub 


AN 


ss 


IC 


AN 


SS 


IC 


AN 


SS 


IC 


1 

238 

0.393 

4.2 

0.8 

0.050 

2 

28,203 

0.278 

3.5 

1.4 

0.049 

3 

2.218,636 

0.227 

3.4 

1.5 

0.050 

1 

179 

0.387 

4.5 

1.1 

0.056 

2 

15,931 

0.274 

4.1 

1.5 

0.056 

3 

939,929 

0.224 

4.0 

1.5 

0.054 

1 

38 

0.288 

5.3 

0 

0.053 

2 

703 

0.204 

4.4 

1.0 

0.054 

3 

8,436 

0.166 

3.9 

1.3 

0.052 

1 

248 

0.283 

3.6 

1.6 

0.052 

2 

30,628 

0.200 

2.9 

2.8 

0.057 

3 

2,511,496 

0.163 

2.6 

3.0 

0.056 

1 

177 

0.278 

4.0 

1.1 

0.051 

2 

15,576 

0.196 

3.7 

2.4 

0.062 

3 

908,600 

0.160 

3.1 

2.7 

0.058 

1 

42 

0.272 

2.4 

2.4 

0.048 

2 

861 

0.192 

2.3 

2.9 

0.053 

3 

11,480 

0.157 

2.2 

2.7 

0.049 

1 

244 

0.351 

0.8 

4.1 

0.049 

2 

29,646 

0.248 

1.3 

3.8 

0.051 

3 

2,391,444 

0.203 

1.5 

3.5 

0.049 

1 

179 

0.392 

0.6 

3.9 

0.045 

2 

15,931 

0.277 

0.8 

3.5 

0.043 

3 

939,929 

0.226 

0.9 

3.8 

0.047 

1 

42 

0.358 

0 

2.4 

0.024 

2 

861 

0.253 

0.9 

3.2 

0.040 

3 

11,480 

0.207 

1.6 

2.8 

0.044 

AN  =  annual;  SS  =  summer  stratification;  IC  =  ice-covered. 


Table  6d.     Comparison  of  Monte  Carlo  simulation  with  the  parametric  confidence  limit.    Number 
of  cases  in  the  simulation  beyond  the  95%  confidence  interval  were  counted  for  Ca. 


Sample 

Number  of 

Monte  Carlo 

Lake 

Season' 

Size 

Cases 

L 

%<(x-L) 

%>(x+L) 

a 

Blue  Chalk 

AN 

1 

201 

0.538 

2.0 

1.5 

0.035 

2 

20,100 

0.380 

2.0 

2.9 

0.049 

3 

1,333,300 

0.310 

2.0 

3.2 

0.052 

ss 

1 

147 

0.560 

2.0 

1.4 

0.034 

2 

10,731 

0.396 

2.2 

3.1 

0.053 

3 

518,665 

0.323 

2.0 

3.3 

0.053 

IC 

1 

34 

0.493 

2.9 

0 

0.029 

2 

561 

0.349 

1.1 

3.4 

0.045 

3 

5,984 

0.285 

1.8 

2.0 

0.038 

Harp 

AN 

1 

196 

0.510 

2.6 

1.5 

0.041 

2 

19,110 

0.360 

1.8 

3.4 

0.052 

3 

1,235,780 

0.294 

1.8 

2.8 

0.046 

SS 

1 

135 

0.484 

2.2 

2.2 

0.044 

2 

9,045 

0.342 

2.1 

2.9 

0.050 

3 

400,995 

0.279 

2.1 

2.6 

0.047 

IC 

1 

33 

0.605 

0 

6.1 

0.061 

2 

528 

0.428 

1.2 

1.8 

0.030 

3 

5,456 

0.350 

1.6 

3.4 

0.050 

Chub 

AN 

1 

199 

0.704 

0.5 

2.0 

0.025 

2 

19,701 

0.498 

1.0 

3.0 

0.040 

3 

1,293,699 

0.407 

1.5 

3.0 

0.045 

SS 

1 

138 

0.698 

0.7 

2.9 

0.036 

2 

9.453 

0.494 

1.0 

3.1 

0.041 

3 

428,536 

0.403 

1.3 

3.2 

0.045 

IC 

1 

38 

0.762 

0 

2,6 

o.ort. 

2 

703 

0.539 

1.7 

1,6 

0.033 

3 

8,436 

0440 

l.b 

1.3 

0.029 

AN  «=  annual;  SS  ■=  summer  stratification;  IC  ■»  ice-covered. 


Table  6e.     Comparison  of  Monte  Carlo  simulation  with  the  parametric  confidence  limit.    Number 
of  cases  in  the  simulation  beyond  the  95%  confidence  interval  were  counted  for  Mg. 


Lake 


Season' 


Sample 
Size 


Number  of 
Cases 


Monte  Carlo 
%<(x-L)      %>(x+L) 


Blue  Chalk 


Harp 


Chub 


AN 


ss 


IC 


AN 


SS 


IC 


AN 


SS 


IC 


1 

200 

0.155 

2.5 

3.5 

0.060 

2 

19,900 

0.110 

3.1 

3.1 

0.062 

3 

1,313,400 

0.090 

3.0 

2.7 

0.057 

1 

146 

0.158 

2.1 

4.1 

0.062 

2 

10,585 

0.112 

2.9 

3.3 

0.062 

3 

508,080 

0.091 

2.7 

2.9 

0.056 

1 

34 

0.131 

2.9 

2.9 

0.059 

2 

561 

0.093 

1.9 

4.2 

0.061 

3 

5,984 

0.076 

2.0 

3.4 

0.055 

I 

200 

0.157 

1.5 

3.5 

0.050 

2 

19,900 

0.111 

2.1 

2.9 

0.050 

3 

1,313,400 

0.091 

2.1 

2.7 

0.048 

1 

135 

0.152 

1.5 

3.0 

0.044 

2 

9,045 

0.108 

1.9 

2.8 

0.047 

3 

400,995 

0.088 

2.0 

2.8 

0.049 

1 

37 

0.163 

0 

2.7 

0.027 

2 

666 

0.115 

2.0 

2.8 

0.048 

3 

7,770 

0.094 

2.0 

3.0 

0.050 

1 

197 

0.146 

2.5 

1.5 

0.041 

2 

19,306 

0.103 

1.8 

2.3 

0.041 

3 

1,254,890 

0.084 

2.1 

2.6 

0.046 

1 

135 

0.152 

2.2 

2.2 

0.044 

2 

9,045 

0.108 

1.7 

2.5 

0.042 

3 

400,995 

0.088 

2.1 

2.6 

0.047 

1 

39 

0.141 

2.6 

0 

0.026 

2 

741 

0.100 

2.9 

0.9 

0.037 

3 

9,139 

0.081 

2.3 

1.4 

0.037 

AN  =  annual;  SS  =  summer  stratification;  IC  =  ice-covered. 


Table  6f.     Comparison  of  Monte  Carlo  simulation  with  the  parametric  confidence  limit.    Number 
of  cases  in  the  simulation  beyond  the  95%  confidence  interval  were  counted  for  Na. 


Sample       Number  of  Monte  Carlo 

Lake  Season'        Size  Cases  L  %<(x-L)      %>(x+L)  a 


Blue  Chalk             AN                1  201  0.236  1.5 

2  20,100  0.167  1.7 

3  1,333,300  0.136  1.7 

SS                   1  148  0.230  2.0 

2  10,878  0.163  1.9 

3  529,396  0.133  1.9 

IC                   1  33  0.226  0 

2  528  0.160  0.4 

3  5,456  0.131  1.2 


Harp                       AN                1  202  0.218  1.0 

2  20,301  0.154  1.6 

3  1,353,400  0.126  1.9 

SS                  1  135  0.195  0.7 

2  9,045  0.138  0.9 

3  400,995  0.113  1.1 

IC                  1  39  0.237  5.1 

2  741  0.167  1.8 

3  9,139  0.137  2.5 


Chub                       AN                1  202  0.243  0 

2  20,301  0.172  0,7 

3  1,353,400  0.141  1.0 

SS                   1  138  0.212  0.7 

2  9.453  0.150  1.5 

3  428.536  0.122  1.7 

IC                   1  41  0.309  0 

2  820  0.219  0 

3  10.660  0.179  0.1 


AN  •=  annual;  S.S  «^  summer  straiificatuin.  IC  -  ice-covered 


2.5 

0.040 

3.2 

0.048 

3.4 

0.051 

2.7 

0.047 

3.7 

0.057 

3.4 

0.052 

3.0 

0.030 

3.8 

0.041 

3.4 

0.045 

3.0 

0.040 

2.9 

0.044 

3.1 

0.050 

4.4 

0.052 

4.3 

0.051 

3.8 

0.049 

5.1 

0.103 

1.8 

0.037 

2.5 

0.050 

5.0 

0.050 

4.1 

0.049 

3.8 

0.048 

4.3 

0.051 

3.3 

0.048 

3.1 

0.047 

7.3 

0.073 

3.9 

0.039 

3.2 

0.033 

Table  6g.     Comparison  of  Monte  Carlo  simulation  with  the  parametric  confidence  limit.    Number 
of  cases  in  the  simulation  beyond  the  95%  confidence  interval  were  counted  for  K. 


Lake 


Season' 


Sample 
Size 

Number  of 
Cases 

L 

1 

2 
3 

201 
20,100 
1,333,300 

0.094 
0.067 
0.054 

1 
2 
3 

149 

11,026 

540,274 

0.098 
0.069 
0.027 

1 
2 
3 

32 
496 

4,960 

0.086 
0.061 
0.050 

1 
2 
3 

201 
20,100 
1,333,300 

0.105 
0.074 
0.061 

1 
2 
3 

135 

9,045 

400,995 

0.097 
0.069 
0.056 

1 
2 
3 

39 

741 
9,139 

0.129 
0.091 
0.074 

1 
2 
3 

202 

20,301 

1,353,400 

0.104 
0.074 
0.060 

1 
2 
3 

138 

9,453 

428,536 

0.103 
0.073 
0.059 

1 
2 
3 

41 

820 
10,660 

0.113 
0.080 
0.065 

Monte  Carlo 
%<(x-L)      %>(x+L) 


Blue  Chalk 


Harp 


Chub 


AN 


ss 


IC 


AN 


SS 


IC 


AN 


SS 


IC 


3.0 
3.1 

2.7 

2.7 
3.1 

2.7 

3.1 
3.6 

3.2 


3.0 
2.3 
2.1 

3.7 
2.4 
2.0 

2.6 
3.5 
3.3 


0.5 
1.6 

1.6 

0.7 
1.2 
1.4 

0 

2.4 

2.0 


2.5 

0.055 

3.3 

0.064 

3.4 

0.061 

3.4 

0.060 

3.9 

0.070 

3.9 

0.066 

0 

0.031 

1.6 

0.052 

2.0 

0.052 

3.5 

0.065 

3.3 

0.056 

3.7 

0.058 

1.5 

0.052 

2.6 

0.050 

3.1 

0.051 

2.6 

0.051 

2.0 

0.055 

1.5 

0.048 

1.5 

0.020 

3.6 

0.052 

2.9 

0.045 

2.2 

0.029 

3.1 

0.043 

2.9 

0.044 

0 

0 

1.5 

0.039 

1.6 

0.036 

AN  =  annual;  SS  =  summer  stratification;  IC  =  ice-covered. 


Table  6h.     Comparison  of  Monte  Carlo  simulation  with  the  parametric  confidence  limit.    Number 
of  cases  in  the  simulation  beyond  the  95%  confidence  interval  were  counted  for  CI. 


Sample 

Number  of 

Monte  Carlo 

Lake 

Season' 

Size 

Cases 

L 

%<(x-L) 

%>(x+L) 

a 

Blue  Chalk 

AN 

1 

204 

0.128 

2.5 

2.0 

0.044 

2 

20,706 

0.091 

2.4 

2.4 

0.048 

3 

1,394,204 

0.074 

2.3 

2.4 

0.048 

ss 

1 

149 

0.118 

3.4 

0.7 

0.040 

2 

11,026 

0.084 

3.2 

1.5 

0.047 

3 

540,274 

0.068 

3.0 

2.1 

0.051 

IC 

1 

36 

0.182 

0 

5.6 

0.056 

2 

630 

0.129 

0 

5.9 

0.059 

3 

7,140 

0.105 

0.4 

4.3 

0.047 

Harp 

AN 

1 

203 

0.243 

3.0 

3.0 

0.059 

2 

20,503 

0.172 

2.7 

2.7 

0.054 

3 

1,373,701 

0.140 

1.9 

2.9 

0.048 

SS 

1 

135 

0.252 

2.2 

3.0 

0.052 

2 

9,045 

0.178 

1.3 

3.1 

0.044 

3 

400,995 

0.145 

1.7 

3.2 

0.050 

IC 

1 

40 

0.243 

2.5 

0 

0.025 

2 

780 

0.172 

4.3 

1.5 

0.057 

3 

9,880 

0.140 

2.8 

1.8 

0.046 

Chub 

AN 

1 

201 

0.239 

1.0 

2.5 

0,035 

2 

20,100 

0.169 

1.7 

2.7 

0.044 

3 

1,333.300 

0.138 

2.0 

2.8 

0.048 

SS 

1 

138 

0.230 

1.5 

2.9 

0.044 

2 

9,453 

0.163 

1.8 

2.8 

0.046 

3 

428,536 

0.403 

1.3 

3.2 

0.045 

IC 

1 

42 

0.298 

0 

24 

0,024 

2 

861 

0.208 

1.2 

28 

0  040 

3 

11.480 

0.170 

1.3 

3.1 

0.045 

AN  >-  annual;  SS  •=  summer  stratification;  IC  -  ice-covered. 


Table  6i.      Comparison  of  Monte  Carlo  simulation  with  the  parametric  confidence  limit.    Number 
of  cases  in  the  simulation  beyond  the  95%  confidence  interval  were  counted  for  SO^. 


Lake 


Season' 


Blue  Chalk 


Harp 


Chub 


AN 


SS 


IC 


AN 


SS 


IC 


AN 


SS 


IC 


Sample 
Size 

Number  of 
Cases 

L 

1 

2 
3 

199 

19,701 
1,293,699 

1.261 
0.892 
0.728 

1 
2 
3 

148 

10,878 

529,396 

1.276 
0.902 

0.737 

1 
2 
3 

32 

496 

4,960 

1.049 
0.741 
0.605 

1 

2 
3 

198 

19,503 

1,274,196 

1.508 
1.066 
0.870 

1 
2 
3 

133 

8,778 

383,306 

1.484 
1.050 
0.857 

1 
2 
3 

38 

703 

8,436 

1.570 
1.110 
0.906 

1 

2 

3 

200 

19,900 

1,313,400 

2.131 
1.507 
1.230 

I 
2 
3 

136 

9,180 

410,040 

2.090 
1.478 
1.207 

1 
2 
3 

42 

861 

11,480 

2.300 
1.626 
1.328 

Monte  Carlo 
5<(x-L)       %>(x+L) 


0.5 
1.9 
2.1 

0.7 
1.9 
2.1 

0 

2.1 
2.5 


1.5 
1.3 
1.6 

2.3 
1.5 
1.7 

0 

1.3 

1.8 


0 

0.5 

1.4 

0 

0.6 

1.4 

0 

0.2 

1.7 


2.0 
2.7 
2.7 

2.0 

2.7 
2.7 

0 
0.1 

2.4 


1.0 
2.8 

2.7 

1.5 
3.1 
2.8 

0 

3.4 

3.0 


2.5 
2.9 
2.9 

2.2 
2.5 
2.9 

0 

3.5 

2.4 


0.025 
0.045 
0.048 

0.027 
0.046 
0.048 

0 

0.022 

0.048 


0.025 
0.041 
0.043 

0.038 
0.046 
0.045 

0 

0.047 

0.048 


0.025 
0.035 
0.042 

0.022 
0.032 
0.043 

0 

0.037 
0.040 


'    AN  =  annual;  SS  =  summer  stratification;  IC  =  ice-covered. 


Table  6j.     Comparison  of  Monte  Carlo  simulation  with  the  parametric  confidence  limit.    Number 
of  cases  in  the  simulation  beyond  the  95%  confidence  interval  were  counted  for  DOC. 


Monte  Carlo 
Lake  Season'        Size  Cases  L  %<(x-L)      %>(x+L) 


Blue  Chalk  AN  1  137  0.445  1.5 

2.0 
1.8 

SS  1  99  0.407  1.0 

1.0 
1.1 

IC  1  23  0.357  4.3 

1.5 
1.9 


Harp  AN  1  144  1.678  1.4 

1.6 

1.3 

SS  1  96  1.739  1.0 

1.7 
1.6 

IC  1  26  0.668  3.8 

2.8 
1.9 


Chub  AN  1  129  0.928  0 

1.0 
1.7 

SS  1  88  0.883  0 

0.5 
1.0 

IC  1  25  0.893  0 

3.5 
2.9 


AN  "-  annual;  .S.S  ■=  suniim'r  strati!  icatinn,  IC  »  ice-cmiMOii 


Sample 
Size 

Number  of 
Cases 

L 

1 
2 
3 

137 
9,316 

491,220 

0.445 
0.315 
0.257 

1 
2 
3 

99 

4,851 
156,849 

0.407 
0.288 
0.235 

1 
2 
3 

23 

253 
11,771 

0.357 
0.252 
0.206 

1 
2 
3 

144 

10,296 

487,344 

1.678 
1.186 
0.969 

1 
2 
3 

96 

4,560 

142,880 

1.739 
1.229 
1.004 

1 

2 
3 

26 

325 
2,600 

0.668 
0.472 
0.385 

1 
2 
3 

129 
8,256 

349,504 

0.928 
0.656 
0.536 

1 
2 
3 

88 

3.828 
109.736 

0.883 
0.624 
0.510 

1 
2 

3 

25 

300 

2.300 

0.893 
0.631 
0.516 

2.9 

0.044 

3.6 

0.056 

3.7 

0.055 

3.0 

0.040 

3.5 

0.045 

3.8 

0.050 

4.3 

0.087 

3.4 

0.049 

1.7 

0.037 

4.2 

0.056 

5.2 

0.068 

6.5 

0.078 

5.2 

0.063 

5.2 

0.069 

6.7 

0.083 

0 

0.038 

1.6 

0.044 

2.6 

0.045 

3.9 

0.039 

2.9 

0.039 

2.9 

0.046 

3.4 

0,034 

3.4 

0.039 

3.3 

0.043 

0 

0 

Ob 

0042 

1.1 

0.040 

Table  6k.     Comparison  of  Monte  Carlo  simulation  with  the  parametric  confidence  limit.    Number 
of  cases  in  the  simulation  beyond  the  95%  confidence  interval  were  counted  for  colour. 


Sample       Number  of  Monte  Carlo 

Lake  Season^       Size  Cases  L  %<(x-L)      %>(x+L) 


Blue  Chalk             AN                1  91  4.636  3.3 

2  4,095  3.278  2.5 

3  121,485  2.677  2.7 

SS                   1  64  15.829  0 

2  2,016  11.192  0 

3  41,664  9.139  0 

IC                   1  17  3.663  0 

2  136  2.590  3.5 

3  680  2.115  2.0 


Harp                         AN                 1  95  8.743  5.3 

2  4,465  6.182  3.2 

3  138,415  5.048  3.0 

SS                   1  65  8.244  6.2 

2  2,080  5.830  3.3 

3  43,680  4.760  3.2 

IC                   1  17  6.771  0 

2  136  4.788  0.3 

3  680  3.909  0.4 


Chub                       AN                1  90  16.910  3.3 

2  4,005  11.957  2.9 

3  117,480  9.763  2.7 

SS                   I  62  19.576  1.6 

2  1,891  13.843  1.6 

3  37,820  11.302  1.6 

IC                   1  18  17.954  0 

2  153  12.696  0.9 

3  816  10.366  2.5 


*    AN  =  annual;  SS  =  summer  stratification;  IC  =  ice-covered. 


2.8 

0.055 

2.2 

0.047 

2.3 

0.049 

1.6 

0.016 

3.1 

0.031 

4.6 

0.046 

5.9 

0.059 

3.1 

0.066 

1.9 

0.038 

3.2 

0.084 

1.9 

0.050 

2.0 

0.049 

0 

0.062 

1.3 

0.045 

1.4 

0.046 

5.9 

0.059 

3.5 

0.038 

3.9 

0.043 

3.3 

0.067 

2.5 

0.054 

2.4 

0.051 

3.2 

0.048 

3.7 

0.053 

4.4 

0.060 

0 

0 

0.9 

0.019 

2.0 

0.046 

Table  61.      Comparison  of  Monte  Carlo  simulation  count  of  number  of  Monte  Carlo  cases  beyond 
the  95%  confidence  limit  (L)  for  total  aluminum. 


Lake 


Season^ 


Sample 
Size 


Number  of 
Cases 


Monte  Carlo 
%<(x-L)      %>(x+L) 


Blue  Chalk 


Harp 


Chub 


AN 


ss 


IC 


AN 


SS 


IC 


AN 


SS 


IC 


94 

4,371 

134,044 

63 
1,953 

39,711 

19 
171 
969 


107 

5,671 

198,485 

68 

2,278 
50,116 

23 

253 
1,771 


111 

6,105 

221,815 

73 

2,628 

62.196 

24 

276 

1.932 


21 
15 
12 

16.33 

11.547 
9.428 

25.434 
17.984 
14.684 


35 

25 
20 

31.833 
22.509 
18.379 

49.657 
35.113 
28.670 


52 
37 
30 

47.584 
33.647 

27.472 

51.42 

36.359 

29.687 


0.9 
0.4 
0.5 

1.5 

0.7 
0.6 

0 
0 
0 


1.8 
2.9 
2.8 

1.4 
2.4 


4.2 
5.0 
4.2 


8.5 

0.085 

6.1 

0.061 

4.6 

0.046 

6.3 

0.063 

4.8 

0.048 

4.3 

0.043 

5.3 

0.053 

4.4 

0.044 

4.6 

0.046 

3.7 

0.047 

4.6 

0.050 

4.9 

0.054 

4.4 

0.059 

5.1 

0.058 

4.8 

0.054 

4.3 

0.043 

8.1 

0.081 

6.4 

0.064 

4.5 

0.063 

4.9 

0.078 

4.0 

0.068 

4.1 

0,055 

4.1 

0065 

3.5 

0062 

0 

0  042 

1.4 

0.064 

1.4 

0.056 

AN  -^  annual;  SS  -  summi-r  stratificaiion,  IC  -  ice-covered 


Table  7.  Estimated  sampling  frequency  required  of  lake  water  chemistry  parameters  for 

various  tolerances  about  the  annual  mean  (x),  with  95%  confidence,  calculated 
using  the  standard  error  of  multi-year  data  sets  for  Blue  Chalk  (BC),  Harp  (HP)  and 
Chub  (CB)  Lakes.   Quality  control  (QC)  replicates  are  from  Locke  (1989). 


Required  Sampli 

ing  Frequency 

Distance 

QC  Replicates  for 

from  X 

Analytical 

DRC  Lakes 

10%  of 

20%  of 

with  one 

sample 

Para- 

Mean 

Mean 

% 

meter 

Lake 

S.D. 

S.D. 

C.V. 

X 

X 

Absolute     ofx 

Cond 

0.22  mS 

0.23  uS 

1% 

BC 

145 

132 

78 

1 

1 

2.65 

9 

HP 

98 

89 

38 

1 

1 

2.18 

6 

CB 

263 

240 

144 

1-2 

1 

3.57 

12 

Alk 

0.05  me/L 

0.114  me/L 

15% 

BC 

209 

40 

2-3 

3-4 

1 

0.72 

18 

HP 

78 

15 

1 

1-2 

1 

0.44 

14 

CB 

69 

13 

16 

35 

9 

0.42 

59 

PH 

0.032  pH 

0.039  dH 

i% 

BC 

151 

102 

36 

1 

1 

0.39 

6 

HP 

78 

53 

20 

1 

1 

0.28 

5 

CB 

120 

81 

38 

1 

1 

0.35 

6 

Ca 

0.151  mg/L 

0.076  mfi/L 

m. 

BC 

12-13 

50 

42 

3-4 

1 

0.54 

20 

HP 

11-12 

45 

30 

2-3 

1 

0.51 

16 

CB 

21-22 

86 

77 

6-7 

1-  2 

0.70 

26 

Mg 

0.042  me/L 

0.018  mg/L 

2% 

BC 

13-14 

74 

109 

4-5 

1-  2 

0.16 

21 

HP 

13-14 

76 

71 

2-3 

1 

0.16 

17 

CB 

11-12 

65 

114 

4-5 

1-  2 

0.15 

21 

Na 

0.07  mg/L 

0.026  mR/L 

3% 

BC 

11-12 

82 

101 

9-10 

2-  3 

0.24 

30 

HP 

9-10 

70 

53 

4-5 

1-  2 

0.22 

22 

CB 

12 

87 

134 

12-13 

2-  3 

0.24 

35 

Table  7.    (Cont'd) 


Required  Sampli 

ing  Frequency 

Distance 

QC  Replicates  for 

from  X 

Analytical 

DRC  Lakes 

10%  of 

20%  of 

with  one 

sample 

Para- 

Mean 

Mean 

% 

meter 

Lake 

S.D. 

S.D. 

C.V. 

X 

X 

Absolute     ofx 

K 

0.04  me/L 

0.01  ms/L 

2% 

BC 

6 

96 

175 

6-7 

1-  2 

0.098 

27 

HP 

6 

96 

82 

3-4 

1 

0.098 

18 

CB 

6 

96 

150 

5-6 

1-  2 

0.098 

25 

CI 

0.051  me/L 

0.045  me/L 

11% 

BC 

6-7 

8 

7 

8-9 

2-3 

0.13 

30 

HP 

20-21 

26 

7-8 

9-10 

2-3 

0.23 

30 

CB 

21-22 

28 

20 

23-24 

6 

0.24 

49 

SO, 

0.167  me/L 

0.141  me/L 

2% 

BC 

28 

40 

43 

1-2 

1 

0.888 

13 

HP 

20 

28 

21 

1 

1 

0.745 

9 

CB 

22 

30 

26 

1-2 

1 

0.776 

10 

DOC 

0.148  mg/L 

0.115  me/L 

4% 

BC 

8-9 

15-16 

39 

6-7 

1-2 

0.45 

25 

HP 

126 

215 

106 

16-17 

4-5 

1.7 

41 

CB 

38 

64 

22 

3-4 

1 

0.92 

19 

Colour 

0.52  units 

2.4  units 

18% 

BC 

80 

3-4 

18-19 

60 

14-15 

4.7 

77 

HP 

283 

13-14 

5-6 

17-18 

4-5 

8.7 

41 

CB 

1,110 

50 

4-5 

13-14 

3-4 

16.9 

36 

Al 

5.3  me/L 

.3,3  mfi/L 

s% 

BC 

16 

41 

519 

333 

83 

21.2 

182 

HP 

32 

84 

76 

49 

i: 

30.0 

70 

CB 

97 

250 

46 

30 

7-  8 

52,1 

54 

Table  8.  Estimated  sampling  frequency  required  of  lake  water  chemistry  parameters  for  various 
tolerances  about  seasonal  (summer  stratified,  SS;  and  ice  cover,  IC)  mean  (x),  with  95% 
confidence,  calculated  using  the  standard  error  of  multi-year  data  sets  for  Blue  Chalk  (BC), 
Harp  (HP)  and  Chub  (CB)  Lakes.    Quality  control  (QC)  replicates  are  from  Locke  1989. 


Sampling  Frequency 


QC  Replicates  for 
DRC  Lakes 


Distance  from  x 
with  one  sample 


Para-  Analytical  Mean  Mean 

meter     Lake     Season  S.D.  S.D.  C.V. 


10%  of  20%  of 


Absolute     ofx 


Cond 


Alk 


PH 


BC 


HP 


CB 


BC 


HP 


CB 


BC 


HP 


CB 


0.22  uS 

0.23  fiS 

1% 

SS 

156 

142 

84 

1            1 

2.74 

9 

IC 

108 

98 

56 

1            1 

2.27 

7 

SS 

46 

42 

17-18 

1            1 

1.49 

4 

IC 

260 

238 

98 

1            1 

3.55 

10 

SS 

240 

220 

133 

1-2          1 

3.41 

12 

IC 

317 
0.05  mR/L 

290 
0.114  mg/L 

164 

15% 

1-2          1 

3.92 

13 

SS 

234 

45 

1-2 

3-4           1 

0.764 

19 

IC 

51 

9-10 

1 

1             1 

0.359 

9 

SS 

68 

13 

1 

1-2           1 

1.49 

4 

IC 

76 

14-15 

1 

1-2           1 

3.55 

10 

SS 

74 

14 

16-17 

38       9-10 

0.431 

62 

IC 

74 
0.032  units 

14 
0.039  units 

16-17 

1% 

37       9-10 

0.429 

61 

SS 

147 

99 

35 

0.388 

6 

IC 

81 

55 

20 

0.288 

5 

SS 

76 

51 

19-20 

0.278 

4 

IC 

72 

49 

19 

0.272 

4 

SS 

124 

83 

39 

0.357 

6 

IC 

126 

85 

41 

0.359 

6 

Table  8.  (Cont'd) 


Sampling  Frequency 


QC  Replicates  for  Distance  from  x 

DRC  Lakes  with  one  sample 


Para-  Analytical  Mean  Mean  10%  of  20%  of  % 

meter     Lake     Season  S.D.  S.D.  C.V.  x  x  Absolute     ofx 


K 

0.04  mE/L 

0.01 

mg/L 

2% 

BC 

ss 

6 

96 

175 

7-8 

1-2 

0.098 

26 

IC 

4-5 

74 

129 

5-6 

1-2 

0.086 

23 

HP 

ss 

4-5 

74 

66 

2-3 

1 

0.086 

16 

IC 

10-11 

167 

138 

5-6 

1 

0.129 

24 

CB 

ss 

7-6 

108 

169 

6-7 

1-2 

0.104 

26 

IC 

8-9 

129 

192 

7-8 

1-2 

0.114 

28 

CI 

0.051  mg/L 

0.045 

mg/L 

11% 

BC 

ss 

5-6 

6-7 

6-7 

7-8 

1-2 

0.118 

28 

IC 

8-9 

11-12 

8-9 

10-11 

1 

0.151 

33 

HP 

ss 

21 

27 

7-8 

9-10 

1 

0.235 

31 

IC 

23 

29 

8-9 

9-10 

1 

0.243 

31 

CB 

ss 

20 

26 

18 

22 

1-2 

0.229 

47 

IC 

25 

33 

21 

25 

1-2 

0.257 

50 

SO, 

0.167  mR/L  0.141 

mg/L 

2% 

BC 

ss 

31 

43 

47 

1-2 

2-3 

0.929 

14 

IC 

19 

27 

28 

1-2 

1-2 

0.733 

11 

HP 

ss 

23 

32 

25 

1 

1 

0.804 

10 

IC 

11-12 

16-17 

12-13 

1 

1-2 

0.572 

7 

CB 

ss 

2! 

30 

26 

1-2 

2-3 

0.770 

10 

IC 

25 

35 

29 

1-2 

4-5 

0.829 

II 

Table  8.  (Cont'd) 


Sampling  Frequency 


Para- 
meter    Lake     Season 


QC  Replicates  for 

Distance  from  x 

Analytical 

DRC  Lakes 

10%  of  20%  of 

with  one  sample 

Mean 

Mean 

% 

eason 

S.D. 

S.D. 

C.V. 

X 

X 

Absolute 

of  X 

0.151  me/L 

0.076  me/L 

3% 
46 

4-5 

1-2 

0.561 

SS 

13-14 

54 

20 

IC 

10-11 

42 

36 

3-4 

1 

0.494 

18 

SS 

10-11 

41 

27 

2-3 

1 

0.484 

16 

IC 

16-17 

64 

40 

3-4 

1 

0.606 

19 

SS 

21 

84 

76 

6-7 

1-2 

0.698 

26 

IC 

25 

0.042  me/L 

101 
0.018  me/L 

85 

2% 

7-8 

1-2 

0.762 

28 

SS 

14-15 

78 

115 

4-5 

1-2 

0.159 

21 

IC 

9-10 

53 

77 

3-4 

1 

0.131 

18 

SS 

13-14 

72 

66 

2-3 

1 

0.153 

16 

IC 

15-16 

86 

77 

3-4 

1 

0.166 

18 

SS 

13-14 

72 

126 

5-6 

1-2 

0.153 

22 

IC 

11-12 

61 

102 

4-5 

1-2 

0.141 

20 

Ca 


Mg 


BC 


HP 


CB 


BC 


HP 


CB 


Na  0.07  me/L  0.026  me/L        3% 

BC 


HP 


CB 


SS 

10-11 

78 

99 

8-9 

2-3 

0.229 

30 

IC 

10-11 

75 

82 

7-8 

1-2 

0.225 

27 

SS 

7-8 

57 

44 

4 

1 

0.196 

20 

IC 

11-12 

83 

56 

5-6 

1-2 

0.237 

22 

SS 

9-10 

66 

108 

9-10 

2-3 

0.212 

31 

IC 

19-20 

142 

185 

16-17 

4-5 

0.310 

41 

Table  8. 


(Cont'd) 


Sampling  1 

Frequency 

Distance 

QC  Replicates  for 

from  X 

Analytical 

DRC  Lakes 

10%  of  20%  of 

with  one 

sample 

Para- 

Mean 

Mean 

% 

meter 

Lake 

Season 

S.D. 

S.D. 

C.V. 

X 

X 

Absolute     ofx 

DOC 

0.115  mg/L  0.115  me/L 

4% 

BC 

SS 

7-8 

12-13 

29 

4-5 

1-2 

0.408 

21 

IC 

5-6 

9-10 

25 

3-4 

1 

0.357 

20 

HP 

SS 

138 

229 

113 

17-18 

4-5 

1.74 

42 

IC 

20 

34 

10-11 

2-3 

1 

0.668 

16 

CB 

SS 

36 

59 

21 

3-4 

1 

0.882 

18 

IC 

36 

60 

18-19 

3-4 

1 

0.894 

17 

Colour 

0.052  H  units 

2.4  H  units 

18% 

BC 

SS 

53 

2-3 

6-7 

31 

7-8 

3.80 

56 

IC 

50 

2-3 

8-9 

101 

25 

3.67 

100 

HP 

SS 

252 

11-12 

7-8 

16-17 

4-5 

8.25 

41 

IC 

169 

7-8 

8-9 

8-9 

2-3 

6.76 

29 

CB 

SS 

912 

43 

18 

12-13 

3-4 
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Fig.  4  Frequency  distribution  of  12  parameters  identified  in  Table  2  in  Blue  Chalk 

Lake.    The  theoretical  normal  distribution  about  the  mean  is  denoted  by 
asterisks. 
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Fig.  5  Frequency  distribution  of  12  parameters  identified  in  Table  2  in  Chub  Lake. 

The  theoretical  normal  distribution  about  the  mean  is  denoted  by  asterisks. 
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Fig.  6  Frequency  distribution  of  12  parameters  identified  in  Table  2  in  Harp  Lake. 

The  theoretical  normal  distribution  about  the  mean  is  denoted  by  asterisks. 
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Fig.  7  The  distance  from  the  mean  value  as  a  function  of  sample  size  (based  on  95% 

confidence  limits)  for  12  parameters  for  each  of  2  seasons  for  Blue  Chalk, 
Chub  and  Harp  Lakes. 
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Fig.  1  Average  concentrations  of  12  key  parameters  identified  in  Table  2  in  Blue 

Chalk  Lake.  Annual  (AN)  and  seasonal  (IF  =  ice-free,  IC  =  ice-covered,  FO 
=  fall  overturn,  SS  =  summer  stratified,  SO  =  spring  overturn)  means  are 
shown. 

a.  Conductivity  e.  Mg  i.  SO,' 

b.  Alkalinity  f.  Na  j.  DOC 

c.  pH  g.  K  k.  Colour 

d.  Ca  h.  CI  1.  Al 

Fig.  2  Average  concentrations  of  12  key  parameters  identified  in  Table  2  in  Chub 

Lake.  Aimual  (AN)  and  seasonal  (IF  =  ice-free,  IC  =  ice-covered,  FO  =  fall 
overturn,  SS  =  summer  stratified,  SO  =  spring  overturn)  means  are  shoN^-n. 

a.  Conductivity  e.  Mg  i.  SO/ 

b.  Alkalinity  f.  Na  j.  DOC 

c.  pH  g.  K  k.  Colour 

d.  Ca  h.  CI  1.  Al 

Fig.  3  Average  concentrations  of  12  key  parameters  identified  in  Table  2  in  Harp 

Lake.  Aimual  (AN)  and  seasonal  (IF  =  ice-free,  IC  =  ice-covered,  FO  -  fall 
overturn,  SS  =  summer  stratified,  SO  =  spring  overturn)  means  are  shown. 

a.  Conductivity  e.  Mg  i.  SO/ 

b.  Alkalinity  f.  Na  j.  DOC 

c.  pH  g.  K  k.  Colour 

d.  Ca  h.  CI  1.  M 
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Fig  .    7.x 
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